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Abstract 


The deleterious effects of fretting on the fatigue properties of a material have been known since 
the work of Warlow-Davies in 1941. However, a widely accepted fretting fatigue life prediction 
method still does not exist and debates persist as to what the critical paramieter(s) for fretting 
fatigue are. This work demonstrates that the surface stresses due to contact can be used to 
characterize the damaging effects of fretting on the fatigue properties of a material. 

The elastic stress analvsis performed by Hamilton for a sphere on flat plane fretting geometry 
is used to determine the stresses due to contact for various experimental conditions provided 
in previous work. It is shown that the specimen fretting fatigue life is directly related to the 
surface stresses induced in the material by fretting. 

Fretting fatigue life prediction methods are analvzed. The stress-life approach. which is 
widely used in industry, is examined for its applicability for fretting fatigue conditions. The 
damage tolerant approach is used to determine the relative magnitudes of the crack nucleation 
and long crack propagation stages of fretting fatigue life. 

An analysis of fretting fatigue palliatives, focusing on shot peening and coatings. is pre- 
sented. A svstematic method to determine the optimum shot peening depth for fretted com- 
ponents is proposed. Tle method uses the elastic stress field expressions derived by Hamilton 
for a sphere on flat surface contact geometry and is demonstrated for Ti-5Al-4V. However, tlıe 
method is general and can be applied for any material and anv contact geometry which allows 
analvtic evaluation of the stress fields. 
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Chapter 1 


Introduction 


1.1 Definitions and Background 


Fretting fatigue generally refers to the degradation of the fatigue properties of a material due 
to cyclic relative displacement between two contacting materials. The most severe impact on 
the fatigue properties occurs when displacement amplitudes are 10-50 jum [21]. The fretting 
action produces considerable wear damage to the material and also produces very large contact 
stresses. The reduction in fatigue life under these circumstances can be quite significant. The 


National Materials Advisory Board [1] gives the following definitions. 


e Fretting . lretting is a surface damage phenomenon occurring on two contacting surfaces 


having oscillatory relative motion of small amplitude. 


e Fretting-Initiated Fatigue . A condition where the material fatigue strength is de- 


graded by the presence of the following: 


1. Surface stress concentrations resulting at sites of fretting pits. 


2. Surface and subsurface stresses resulting from rubbing friction contact combined 


with the stresses resulting from the overall fluctuating stress field. 


There are many practical examples of structures in which fretting occurs. ‘These structures 
include: bolted and riveted joints, suspension cables, metallic coil ropes, kev-way shaft cou- 


plings. coil wedges of generator rotors and blade dovetail contact sections of turbine engines 


10 


[2]. The fretting problem is made worse by the fact that many of the structures that undergo 
fretting cannot be easily inspected. The first detection of a crack mav be when the structure 
fails. 

The first reported observation of the occurrence of fretting during fatigue testing was made 
bv Eden. Rose, and Cunningham [3] in 1911. They observed corrosion in the holders of the 
fatigue test specimens. The corrosion was such that they had difhculty removing the broken 
test specimens from the testing apparatus for the heavily loaded tests which lasted for several 
days. Eden, Rose. and Cunningham [3] also recorded the first attempt of the use of palliatives 
for fretting fatigue. They oiled the fatigue specimens before putting them in place, however 
they reported that this was only partially effective. The specimens were easier to remove. but 
corrosion products were still produced. The tests conducted by Eden et al [3] were not fretting 
fatigue tests, the observed fretting was noted but not analyzed. 

According to Hoeppner [4] the first systematic investigation of fretting was performed by 
Tomlinson in 1927, in which he determined that surface damage was related to the relative 
displacement between the contacting surfaces. Warlow-Davies [5] investigated the impact of 
fretting on the fatigue properties of a material in 1941 and demonstrated that fretting reduced 
fatigue strength. He did not perform fretting and fatigue tests simultaneously, but fretted the 
specimens and then performed fatigue tests on the fretted specimens. Warlow- Davies performed 
the tests on medium-carbon steel and for a nickel-chromium-molvbdenum alloy steel. The loss 
in fatigue strength after fretting was found to be 13 percent and 18 percent respectively. By 
separating the fretting and fatigue portions of the test, Warlow-Davies was able to measure 
the reduction m fatigue strength due only to the surface damage caused by fretting. Later 
research, such as that conducted by Waterhouse [29], showed that simultaneous fretting and 
fatigue resulted in much higher reductions in fatigue strength. There have been numerous 
other experiments demonstrating the deleterious effects of fretting on fatigue strength. Birch 
[11], Szohwinski & Farris [30] and Adibnazart & Hoeppner [15] are just a few of the other works 
that have also demonstrated the damaging effects of fretting on fatigue life. 

Much work has been performed in fretting fatigue since the investigation conducted by 
Warlow-Davies. This thesis will not attempt to reference all the major work performed in 


the area of fretting fatigue. Hoeppner in [4] provides a historical review on fretting fatigue 





specimen 





Fretting pads 


Figure 1-1: Simplified schematic of bridge type fretting pads. 


research as well as a brief review of the current state of knowledge, which 1s part of an entire 
work devoted to the topic of fretting fatigue. Birch [11] also provides a brief summary of 


fretting research. 


1.2 Motivation and Objectives 


Although the damaging effects of fretting on the fatigue properties of a material have been 
reported quite extensively since the work of Warlow-Davies. there is still not a consensus on 
the critical parameter(s) for fretting fatigue. Waterhouse [9] states that the main factor in 
producing fatigue failure is the additional shear stress in the contact region caused by the 
frictional forces between the contacting bodies, while Nishioka and Hirakawa [10] stress the 
importance of the relative slip amplitude and Adibnazari and Hoeppner [15] demonstrate the 
importance of the normal pressure on fretting fatigue life. According to Birch [11], Beard lists 
more than fiftv parameters of concern for fretting fatigue. The lack of consensus on the critical 
fretting parameter(s) is likely due to the experimental configurations that have been used in 
the past. Many of the past fretting fatigue experiments have been performed using bridge type 
pads as shown in Figure 1-1. 

The bridge pads produce stress singularities at the contact edges making a stress analysis 
very difficult. Many of the experimental configurations in the past also lacked the ability 


to measure all of the fretting parameters, specifically the oscillating tangential contact load. 





which is a verv important fretting parameter. If the oscillating tangential load is not present. 
the problem becomes an indentation problem and not fretting. This inabihty to measure the 
contact load conditions again limited the information available for analysis. The first objective 
of this work is to determine the critical fretting fatigue parameter(s). 

The limitations discussed above are not in this work. The work of Birch [11] and the ongoing 
experimental work of Conner [12] provide the experimental data used in this work. The testing 
apparatus used by Birch and Conner employed spherical fretting pads. The spherical fretting 
pads allow analytical expressions for the stress field due to contact to be determined. which 
was first done by Hamilton and Goodman [13]. ^ The testing apparatus used by Birch and 
Conner. described in detail in [25]. allows measurement of all the fretting parameters required 
to determine the stress fields due to contact. The expressions given in [13] were expressed again 
later bv Hamilton [14] in a more convenient explicit form. The form developed by Hamilton 
[14] is implemented in this work. Hamilton and Goodman provide the elastic stress fields. 
This does not present a problem in this work as the stresses remained in the elastic region for 
all but a few tests. 

The expressions given by Hamilton [14] are for a global sliding contact, which is not the 
condition for most of the experiments analyzed in this work. Most of the tests analvzed in 
this thesis have shp/stick conditions. Chivers and Gordelier [16] show how the Hamilton 
expressions must be modified to accommodate slip/stick fretting conditions. This modification 
has been implemented in this work for the shp/stick experimental conditions when necessary. 
It will be demonstrated in this work that the critical fretting parameters can be captured by 
one parameter, the maximum von Mises effective stress at the surface due to contact. 

The second objective of this work is to analyze fretting fatigue life prediction methods 
using past experimental work. A widely accepted fretting fatigue life prediction method does 
not yet exist. The stress life approach. which is widely used in industry today, is examined 
for its applicabilitv for fretted components. There have been several life prediction methods 
proposed for fretting fatigue, but none have utilized the stress-life approach. Edwards [6] in 
1981 proposed a fretting fatigue life prediction method using linear elastic fracture mechanics. 
Faanes and Fernando [7] in 1994 also proposed a life prediction method using linear elastic 


fracture mechanics, but added a short crack correction to account for the effect of the high 





stress gradient due to fretting on short crack growth. The success of linear clastic fracture 
mechanics applied to fretting has been limited by the difficulty in developing stress intensity 
factors to be used to predict short crack growth rates. Giannakopoulos, Lindley and Suresh 
[46] propose a method which circumvents this problem for sharp edged contacts in their crack 
analogue model. Szolwinski and Farris [30] in a recent paper propose a fretting fatigue life 
prediction method that uses the Smith-Watson-Topper expression [8] to predict the number of 
cycles to crack nucleation and employs linear elastic fracture mechanics to predict the long crack 
growth rate, thereby avoiding the short crack growth rate problem. This work will examine 
both the stress-life approach and the damage tolerant approach used to predict fatigue life. 

The final portion of this work analvzes fretting fatigue palliatives; coatings and shot peening. 
Shot peening has been shown by Waterhouse [38], Bignonnet [39], and Chivers and Gordelier 
[41] to be an effective fretting fatigue palliative. However, a detailed analysis of the shot peening 
process applied to fretted components with the objective of determining the optimum peening 
depth has not vet been performed. [retting produces large tensile and compressive stresses at 
the contact edges and both must be considered when determining the optimum peening depth. 
A systematic method to determine the optimum peening depth for a fretted component is given 
in this work, using the elastic stress field expressions provided by Hamilton. ‘This analvsis also 
includes the effects of changing the coefficient of friction on the optimal peening depth. The 
elastic stress fields due to contact will also be used to provide a brief discussion of tlie benefits 
of protective coatings for fretting fatigue. 

This thesis is organized in the followmg manner: Chapter 2 provides a description of the 
sphere on flat plane contact mechanics used in this work. Chapter 3 uses the works of Hamilton 
[1-4]. and Chivers and Gordelier [16] to analyze the experimental results provided by Birch [11], 
and Conner [12]. Chapter 4 analyzes fretting fatigue life prediction. Chapter 5 discusses 
fretting fatigue palliatives and provides a method to determine the optimum peening depth for 
a fretted component. Chapter 6 provides the conclusions from this work and suggestions for 


further study in work related to fretting fatigue. 





Chapter 2 


Review of Sphere on Flat Plane 


Contact Mechanics 


2.1 Introduction 


Before the analysis of Hamilton [14] and Chivers and Gordelier [16] is applied to experimental 
tests, a brief review of the sphere on flat plane contact mechanics is presented. ‘The expressions 
reviewed here serve as the foundation for the analysis performed in this work. The discussion 


begins with the analysis conducted by Hertz in 1882. 


2.2 Hertzian Contact 


Over one hundred years ago Hertz analyzed the problem of normal, frictionless contacting 
spheres. Hertz solved for the contact pressure distribution, p(r) and the circular radius of 
contact, a, by solving the integral equation governing the boundary conditions. The radius of 


contact as determined by Hertz [17] is given by: 


BP RAR E an | 
qc MERE: 2 (23) 
A( 12, + lt) | TE) 7 E» 





where E and v are Young’ s Modulus and Poisson's ratio respectivelv and. Ris the radius of 


curvature. The indices 1,2 refer to the two contacting bodies. The solution derived by Hertz 
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mE 


(b) 


Figure 2-1: (a) Fretting contact between a shpere and a flat surface. (b) Top view showing 
circular contact region and stick zone, c. 


for normal contacting spheres can be applied to a sphere on a flat plane geometry by noting 


that gl = gi t gl, Rpıane = 00, therefore R = Reyy = Rephere- Huber [18] in 1904, 


Rsphere Rolane ! 


was the first to solve for the stress fields beneath the contacting surfaces. 


2.2.1 Sliding Elastic Sphere on a Flat Surface 


The indentation problem solved by Hertz and Huber becomes a fretting problem when an 
oscillatory tangential force is apphed to the sphere. ‘The problem geometry is shown in Figure 
2-1. 


Once the tangential force, Q. is introduced, sliding or partial slip is also introduced. The 


Amonton-Coulomb law relates the normal pressure and the tangential tractions, q bv: 


lal _ 101 _ 


Where jz is the coefficient of kinetic friction. 


As given first by Mindlin [19] in 1949 and again later by Hamilton [14] the stresses at the 





contact surface for incipient sliding conditions are given by: 





—3P ; e | 
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ut 


Zus and the boundary stresses outside the circular area of contact are zero. As mentioned 
above, these boundary conditions only apply under conditions of gross sliding, that is when 
Q = uP. The expression given by Equation 2.3 is also the expression for contact pressure due 
to indentation without sliding and the expression given for the contact radius by Equation 2.1 
is also unchanged by the addition of a tangential load. Therefore, the presence of the shear 


stresses on the surface has not affected either the circular contact radius or the contact pressure. 


2.2.2 Partial Sliding 


When the apphed tangential force is less than the product of uP, macroscopic sliding does not 
occur, however frictional traction at the contact interface leads to small relative displacement 
over part of the interface. This small relative motion is commonly referred to as micro-slip. 
The remaining portion of the interface deforms without anv relative motion. ‘The region in 
wluch no relative motion occurs is referred to as the stick region. All tlie points on the surface 
within the stick region undergo the same tangential displacements. The stick region is shown 
in Figure 2-1. Mindlin [20] has shown that a stick region exists in the contact interface and 


the radius, c, of the outer boundary of the stick region to be given by: 
C | S 
- = |1 - — Co (2.5) 


As shown in Figure 2-1. the stick region is concentric with the contact area, this geometric 
relationship is also observed for experimental specimens. 
Just as the boundary conditions had to be reestablished for sliding, they must again be 


modified to account for micro-slip. The stresses at the contact surface are now given by [19]: 
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The expression for the normal pressure is the same as in the global sliding condition. The 


relative taugential displacements are given by [19]: 
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2.3  Hamilton's Analysis 


As mentioned in the previous section, the stress fields beneath the surface of a circular indenter 
were first derived by Huber in 1904. However. the elastic stress fields beneath the surface of a 
circular sliding contact on a flat surface were not derived until 1966 by Hamilton and Goodman 
[13]. Mindlin's work in circular sliding contact focused primarily on surface compliance and 
tle micro-slip problem. He did not analyze the stresses below the contact surface. 

It must be stressed that the equations given by Hamilton are for a global shding sphere 
and not for stick/shp conditions. It was shown earlier that the boundary stresses are not 
the same for the two conditions. The stress equations for the global sliding condition must 
be modified slightly to obtain the correct stresses for situations involving micro-slip. Chivers 
and Gordelier [16] demonstrate how the expressions must be modified. The results of a stress 
analvsis performed using the expressions given by Hamilton are shown in Figure 2-2. The 
expressions used to generate this plot are given in Appendix D. The material used to generate 
this plot was Al 7075- T6. Some material properties for tliis alloy are listed in Table 2.1 [11] 
[21] [22]. Spherical fretting pads with a one inch radius were used and global sliding is assumed. 
Using the values for 44, E, and v in Table 2.1 and the elastic stress field equations given by 
Hamilton, the stresses in the x and v directions have been plotted in Figure 2-2 along with 
the von Mises effective stress, Ceg, at z = 0, y — 0. "The stresses have been normalized by 
the maximum contact pressure, [. [ = AE. The von Mises effective stress parameter, c.g. 
has been used because this is a multiaxial problem. The effective stress is used to convert the 


multiaxial stress state to a single stress parameter. The expression for the effective stress is 





0 a 
0.3 

469 MPa 
1.2 

579 MPa 
179 MPa 
1317 MPa 


Material Property 


Young's modulus (E) 
Poisson ratio (1) 
Monotonic yield stress (0 y) 
Coefficient of kinetic friction (p) 
Tensile strength (ors) 
Endurance limit (Cena) 

Fatigue strength coefficient (c) 














































l'atigue strength exponent (b) -0.126 
Fatigue limit stress (op) 230 MPa 
Fatigue ductility coefficient (&) | 0.19 
Fatigue ductility exponent (c) -0.56 


Table 2.1: Properties of 7075-T6 aluminum alloy 
given by: 


1/2 
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na a meee aoo (2.9) 


The stress fields displayed in Figure 2-2 are not the only non-zero stress fields at the surface; 
the shiear stress 7;, and the normal stresses in the z-direction are also non-zero, however they 
both go to zero at the contact edges. The stresses below the surface at a depth of 20 jum 1s 
plotted in Figure 2-3. Comparing Figures 2-2 and 2-3, it can be seen that the stresses due to 
contact quickly decay away, especially the tensile stress at the trailing edge of contact, r = —a, 
in the x-direction. From Figure 2-2 it can be seen that there is a large peak in cy at the trailing 
edge of contact. This stress is the highest stress on the plot. In Figure 2-3, the peak tensile 
stress at the trailing edge is much lower in magnitude and the slope of the a, curve is also much 
smoother. 

The maximum von Mises effective stress is also at the surface. ‘This is true for a sphere 
sliding on a flat surface when u > 0.3. When 4 — 0, or when the tangential force. Q. is zero. 
the maximum von Mises stress shifts below the surface to the point z = 0.5a, for v — 0.3 [14]. 
However, for a sphere fretting on a flat plane. the maximum von Mises and the maximum tensile 
stresses generallv occur at the surface. Although the maximum tensile stress has a severe peak 
at the trailing edge of contact at the surface, the von Mises effective stress 1s roughly constant 


over a verv large region as can be seen in Figure 2-2. 
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Figure 2-2: Contact stresses at the surface due to a sphere sliding on a flat surface for Al 
1015- T6. 
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Figure 2-3: Contact stresses due to a sphere sliding on a flat surface for Al 7075-T6. z = 20 
jum. 
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The analvsis performed in this work will focus on tlıe surface stresses, at the trailing edge of 
contact, with the exception of the shot peening analysis, which analvzes below the surface and 
inside the contact edges. The contact stresses are maximum at the surface and the maximum 
tensile stress occurs in the material at the trailing edge of the contact; therefore. this work 
will focus its analysis on this location. The von Mises effective stress is just as high at other 
locations on the surface, however this is in general due to large compressive stresses at that 
location. It is assumed that the large tensile stress is responsible for crack nucleation. The 
contact stresses given by Hamilton can be greatly simplified for r = —a, z =U and y = 0. the 
problem becomes a biaxial one. The only non-zero stresses at this location are 0, and oy, and 


are given by: 
(1—2v) 4+v 








oo a $ and (2.10) 
214 — 1 37 
gc ale ccu )2 11 


3 8 
Note carefully the position of the signs. The maximum tensile stress in both the x and v 
directions occurs at the trailing edge of contact. Again, it must be emphasized that the 
stresses given in equations 2.10 and 2.11 are for full sliding and not for stick/slip conditions. 
Figure 2-4 displays the curves for e, under incipient sliding and stick/slip conditions. As can 
be seen in this figure. when stick/slip is accounted for, o, is reduced for the given contact load 
conditions. This figure was generated using a c/a ratio of 0.5. 

Recall that the ratio of c/a indicates how close to full sliding or full stick the fretting 
condition is. As c/a — 0, global sliding is introduced, when c/a = 1, full stick conditions 
exist. From Figure 2-4. it can be seen that the magnitude of the maximum tensile stress 
is reduced when micro-slip is taken into consideration. This fact will become more relevant 
later during the shot-peening analvsis of this work. The reduction in stress under stick /slip 
conditions is in direct conflict with experimental results. Experiments have shown that fretting 
is most damaging to the fatigue properties of a system under micro-slip conditions with relative 
displacements of 10-50 ym. The introduction of global sliding improves the fretting fatigue 
life. It has been postulated that this improvement under sliding conditions may be due to 
the wearing away of fatigue cracks from the surface, so that they do not become propagating 


fatigue cracks. This is one example of the limitations of the stress analysis. 
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l'igure 2-4: Plot demonstrating the effects of stick-slip conditions on oz for A} 7075-T6 (c/a = 
0.5). 


In most practical structures and scientific experiments, the tangential load, Q, which changes 
the problem from one of indentation to fretting. is caused by an applied cyclic bulk stress. 
Therefore, to perform a complete analvsis, the stresses induced in the material by the bulk load 
must also be accounted for. If the material behaves elastically, the principle of superposition 


can be applied. Using superposition, the stresses at the surface, at r = a, are given bv: 


Oxtot = Oxc + Oxb and (2 12) 


gl C (213) 


where Croit and Geer are the total stresses in their respective directions. cy, and Oy, are 
the contact stresses in their respective directions andio e Gy» are the bulk stresses in the x 
and v directions respectively. The condition described above, that of a cvclic bulk load in 
combination with the apphed contact loads, was the method used to produce the experimental 
results analyzed in this work. The bulk load also produces eccentricity in the contact geometry. 


The contact area and stick zone are no longer perfectly circular. This effect is discussed in 


215 





detail in [23]. 
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Chapter 3 


Experimental Results 


3.1 Introduction 


As shown earlier, the elastic stress fields for a sphere sliding on a flat plane are known and can 
be found as functions of the normal and tangential loads. ‘The testing apparatus developed at 
the Massachusetts Institute of Technology gives the user the ability to control and measure both 
the normal and tangential contact loads. This apparatus has been used to conduct fretting 
fatigue tests on Al 7075-T6 and Ti-6Al1-4V. The testing apparatus is described in [25]. A 
simplified diagram of the test setup 1s shown in Figure 3-1. These tests provide the necessary 
information to calculate the elastic stress fields due to contact, which will be discussed in detail 
in the following sections. The elastic stress fields at the surface due to contact do not completelv 
describe the fretting phenomenon, there are some shortcomings. ‘The stress analvsis does not 
account for the condition of the surface and may not give the correct location of crack initiation. 
However, the stress analvsis provides a very good starting point for a systematic analysis of the 


problem. 


3.2 Al 7075-T6 Experimental Results 


The experimental data provided in the work performed by Birch [11] and Conner [12] is shown 
in Table 3.1. Table 3.1 provides all the necessary information to calculate the stress fields. The 


bulk stress, ep, for these experiments was applied in the x-direction only, with a stress ratio, 
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Applied Bulk Load 


Figure 3-1: Simplified diagram of the fretting fatigue testing apparatus. 


Care must be taken when superposing the bulk stress. As discussed earlier the trailing 
edge of contact is the point of interest for this analysis; however, during each complete stress 
cycle the trailing edge becomes the leading edge as the bulk stress is reversed. The trailing 
edge during the tensile portion of the bulk stress loading cycle becomes the leading edge during 
the compressive portion of the loading cycle. This location on the specimen undergoes the 
highest tensile stresses and the largest stress amplitude. The opposite edge of contact endures 
lower stresses because it becomes the trailing edge of contact during the compressive portion 
of the bulk loading cvcle. It would be expected that the cracks would initiate from the edge 
of contact which undergoes the highest stresses. This fact is verified in [23] which displays a 
fretting scar with cracking visible only on one side of the circular contact region. 

The loading conditions displayed in Table 3.1 were used to calculate the surface stresses 
at the trailing edge of contact for each test specimen. To better analyze the results of the 
stress analysis, this work has isolated two series of experiments, which correlate to two different 
loading conditions. The first series of tests includes those tests where the bulk stress was 


roughly constant (c483MPa), while the contact loads were varied. Therefore, the effects of the 


OCH Nr (Cycles) 


449.5500 
480,000 
395.000 
361,000 
001,000 
530,000 

3,300,000 (no failure) 
10.680.000 (no failure) 
803.000 
2,940,000 
3,450,000 (no failure) 
616,000 
1,410,000 (no failure) 
2,186,000 (no failure) 
2,500,000 (no failure) 
2,666,000 (no failure) 
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2.005.000 (no failure 
2.608.000 (no failure 
549.000 
516,000 
1,540,000 (no failure) 
2,940,000 (no failure) 
1,777,000 (no failure) 
297.000 
135,000 
480,000 
919.000 
3,050,000 
390,000 
258,063 
181.118 
2,497,226 
3,200,845 (no failure) 
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Table 3.1: Al 7075-T6 fretting fatigue experiments 
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0.82 248 10.680.000 (no failure) 
0.69 264 3.300.000 (no failure) 
0.63 264 3.450,000 (no failure) 
[S 266 2,940.000 
I 280 2.608.000 (no failure) 
0.58 285 616.500 
0.14 202 803.000 
0.73 309 530.000 
0.83 317 480.000 
0.74 334 449,500 
0.72 dod 549.000 
0.58 391 395,000 


0.0 366 361.000 


Table 3.2: Al 7075- T6 Constant Bulk load data analysis 


bulk stress are the same in each test and the difference in life is due only to the difference in 
contact stress. The second series of tests analvzed are made up of tests with constant contact 
loading (P = 20 N, Q = 15 N). while the bulk stress is varied. which isolated the effect of the 
bulk stress on fatigue life. 

The results of the stress calculations for the constant o, series of tests are shown in Table 
3.2. Table 3.2 displavs the calculated values for the contact radius, a, the ratio of the stick 
zone radius to the contact zone radius (c/a), the maximum tensile stress at the trailing edge of 
contact, Ortor and the effective stress due to contact at the trailing edge. Cec. 

Table 3.2 has been arranged in order from lowest to highest a... Note that the fatigue 
lifetimes as a result are nearly arranged in order from greatest number of cycles to lowest. 
This fact indicates that a relationship exists between the surface stresses at the trailing edge of 
contact and the fretting fatigue life of the specimen. This relationship will be discussed further 
in the next chapter, where fretting fatigue life prediction methods are discussed. The fatigue 
life versus effective stress due to contact is plotted in Figure 3-2. Figure 3-2 again demonstrates 
the existence of a relationship between the fatigue life and tlie surface stresses due to contact. 
Recall that each one of these tests were conducted with approximately the same applied bulk 
load. therefore the difference in fatigue life is due to the difference in contact loading only. It 
can be seen that as the contact stresses are reduced, the fatigue life increases. For this bulk 


loading condition it appears that specimen failure will not occur if Gee is maintained less than 
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Figure 3-2: Plot of effective stress due to contact versus number of cvcles to failure for o; — 53 
MPa. 


about 264 MPa, which provides an endurance type stress combination. Also note that in 
none of the tests was the effective stress due to contact greater than the yield strength of the 
material; however, in test four, tlıe total stress in the x-direction does slightlv exceed the vield 
strength of the material. 

The stress values are not tlıe only important values displayed ın Table 3.2; the values 
calculated for a and c/a are also quite important. Thev provide analytical values that can be 
easily checked against the actual scars on the specimens, which is done later in this chapter. 
This gives insight into the validity of the tests. 

The value of c/a not only provides a means to verify testing procedures, it also is a key 
parameter in characterizing the nature of the fretting condition. Recall that as c/a goes to 
zero, global sliding is introduced and when c/a = 1, complete stick is occurring. It can be 
seen from Table 3.2, that in Test 4 global sliding should be occurring. The c/a value can also 
provide insight into the existence of a normal pressure threshold discussed by Adibnazarı and 
Hoeppner [15]. Table 3.3 displays tests with the same value of tangential load (15N), with 
increasing values of normal load. It can be seen that as the normal load is increased, the 


maximum tensile stress and the effective stress due to contact actually decrease. Therefore. for 








Table 3.3: Stress variation as the normal load changes for a constant maximum tangential load. 


a given tangential load. as tle normal load is increased, the contact stresses at the trailing edge 
of contact decrease, and from Table 3.2 and Figure 3-2 the fretting fatigue life should increase. 
This provides a rational explanation of wliv a pressure threshold exists. ‘The stresses decrease 
as the normal load is increased due to the lowering of c/a. As c/a goes to zero, the relative 
displacements also go to zero and complete stick results. 

The results for the second series of tests, in which the contact loads were maintained constant 
while the bulk load was varied. are shown in Table 3.4. The data in Table 3.4 is presented 
in order of ascending bulk load. The values for contact radius, a and c/a are unchanged for 
each test, since the contact loads are identical for each test. The values of a and c/a were 
found to be 213 um and 0.72 respectively. The value of c/a was <1, therefore all the tests 
were conducted under stick/slip conditions. The most interesting aspect of this set of tests 
is the appearance of a fretting fatigue endurance type limit for the given contact conditions. 
It appears that when the bulk stress is reduced below approximately 70 MPa for the given 
contact loading conditions, specimen failure will not result. Figure 3-3 is a plot of the bulk 


tress versus the number of failures to cycle (S-N curve) for these contact conditions. This 


fi 


is not a curve with universal applications. It onlv applies to the contact conditions used in 
this set of tests. A different set of contact load conditions will result in a different endurance 
strength. However, combining the data from the first set of Al 7075-T6 tests with this set 
of data and the known endurance strength of the material a curve that does have universal 
applications can be proposed. This curve is shown ın Figure 3-4. In this figure the so called 
knock-down factor, which is the factor reduction of the endurance strength of the material due 


to fretting, is plotted versus the effective stress due to contact at the trailing edge. 


The curve displaved in Figure 3-4 is a very rough approximation of the reduction in the 
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Table 3.4: Al 7075-T6 stress analvsis for constant contact loading 
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1,540,000 (no failure) 
1.711.000 (no failure) 
2,940,000 (no failure) 
3.200.845 (no failure) 
1.410.000 (no failure) 
516,000 
297.20 
3,050,000 
549,000 
390,000 
258,063 
181,118 
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Figure 3-3: Number of cycles to failure versus op, S — N curve for constant contact load 
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Table 3.5: Data used to plot the knock-down factor curve. 


endurance strength of the material due to fretting. The knock-down factor is given by: 


d Cend nO fretting 
Kr end! g) 
o end (fretting) 


where Cena is the endurance strength of the material. Figure 3-4 was constructed using only 
three data points; several more series of tests at different bulk loads need to be performed to 
determine the accuracy of this plot. However, the accuracy of this plot is not the critical issue, 
the premise on which the plot was constructed is much more important. The plot indicates 
that AF is not a single number for a given material, but it can be determined if the contact 
stresses are known. The fact that the strength reduction changes for various contact loads can 
easily be seen in Figure 3-2. This idea will be examined closer in the next chapter. The data 
points used to make Figure 3-4 are displaved in Table 3.5. When the contact stresses go to 


zero, the knock-down factor goes to one. 


3.3 Ti-GAl-4V Experiments 


The second material analyzed in this work is Ti-6A1-4V [12]. The same test apparatus discussed 
earlier was also used in the Ti-6A1-4V fretting fatigue tests, utilizing the same sphere on flat 
surface fretting geometrv. Table 3.6 provides the necessary material properties to perform the 
stress calculations for this material. Fuchs [22] provides fatigue properties for many different 
materials and was used along with [11] and [45] to produce Table 3.6. 

The loading conditions used to conduct the tests on the titaniuin allov are presented in 
Table 3.7. The results of the stress analysis are displayed in Table 3.8. The value used for the 
coefficient of friction was 44 — 0.5 for the displayed results. The fretting pad radius. R, used in 
these tests was 0.5 inches versus one inch used for the A] 7075- T6 tests. 


A plot of the fatigue life versus the effective stress due to contact is shown in Figure 3-5. 
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Figure 3-4: Plot of the knock-down factor versus effective stress due to contact for Al T075- T6. 


Material Property 


Young's Modulus (E) Mor Gra 
Poisson ratio (1) | 
Coefficient of kinetic friction 0.4-0.5 
Monotonic yield stress (a, ) 925 MPa 
Endurance strength (ena) 525 MPa 
Fatigue strength coefficient (cj) | 1933 MPa 
l'atigue stength exponent (b) -0.10 





Table 3.6: Ti-6Al-4V material properties. 


691.612 
109,574 
107,879 
31.540 (interrupted) 


304,456 
474,656 
143,683 
258. 316 
3.251.000 (no failure) 





Table 3.7: Test parameters for Ti-6A]-4V fretting fatigue tests 


RD 





Sua (MPa) | 0. (MPa) 
1 0 
0.39 
Dio 
0.78 
0.80 


0.75 





Table 3.8: Ti-6Al-4V stress analysis 


This figure includes tests 5-9, where the bulk stress was = 300 MPa. Again, the fatigue life 
increases as the contact stresses are reduced. An endurance type limit is also seen for this 
material. For o, ~ 300 MPa. it appears specimen failure will not occur if o.. Is maintained 
less than about 390 MPa. More experimental tests are needed to confirm this observation. 
Not enough experimental data was available to plot a grapli of the fretting fatigue life versus 


apphed bulk stress. 


3.4 Experimental Observations 


Observations of the fretting pads and specimens have been made using a scanning laser micro- 
scope and an optical microscope. The pictures provided m this section were produced using 
Images provided by the scanning laser microscope. The scanning laser microscope made it 
very simple and quick to perform accurate measurements of the fretting parameters, such as 
the contact radius. a and the radius of the stick zone, c. 

The expressions for the contact radius. a. and the stick zone radius, c, as derived by Hertz 
[17] and Mindlin [19] respectively are given in Chapter 2 of this work. The analytical values 
for these parameters have been calculated for each test conducted on Al 7075-T6 and Ti-6Al- 
4V and will be compared with the actual specimen values when possible. The results of the 
comparison are given in the next section. 

A second objective of this observational work was to determine the location of crack initia- 


tion. To achieve this objective the spherical fretting pads and the unbroken fatigue specimens 
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Figure 3-5: Plot of effective stress due to contact versus number of cycles to failure for Ti-6A1-4V 
is 300 NIPa). 


have been analyzed. The fretting pads provide an opportunity to see the influence of the con- 
tact stresses alone on crack initiation. due to the fact that they are not subjected to an applied 


bulk stress. 


3.4.1 Al 7075-T6 


Not all of the specimens from the 33 Al 7075-T6 tests shown in ‘Table 3.1 were available to be 
analvzed, only about one-half of the tests had pads or specimens that could be analyzed. The 
fretting pads provided nearly all of the observations because when the specimens failed they 
usually failed at the fretting scar. It was very dificult to determine the fretting parameters a, 
and c/a for many of the pads and in some cases it could not be done. The stick and slip zones 
could not be distinguished for manv of the pads that correspond to specimen failure. ‘Table 
3.9 shows the comparison between the predicted (Pred. in the table) and actual (Act. in the 
table) values for the fretting geometry, for those specimens where the comparison was possible. 
L and R correspond to the left and right fretting pads respectively. Test 33. which did not 


ae 


produce specimen failure also displays the values for the specimen and are denoted by “(spec) 
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Figure 3-6: Picture of the left fretting pad used for test 33 with the approximate size of the 
stick zone shown (Al 7075-T6). 


in the table. It can be seen from this table that the actual contact radius is in general larger 
than that predicted by Equation 2.1. In some cases the difference is quite substantial. The 
actual c/a ratios are fairly close to the predicted values, indicating that the value used for u is 
accurate. 

As discussed earlier, a scanning laser microscope was used to observe these specimens and 
to produce images of the scars. Figure 3-6 shows the fretting scar on the left pad used for test 
33. The stick and slip regions are clearly evident and the eccentricity caused by the apphed 
bulk can also be seen. The corresponding fretting fatigue specimen is shown in Figure 3-7. The 
fretting scar on the specimen is not as visible as the scar on the pad. ‘The fretting geometry 1s 
not as easy to distinguish for specimens that were cycled to failure. It is dificult to determine 
even the contact radius of the pads for the failed specimens and nearly impossible to determine 
the stick zone radius. 

The other objective of this section was to determine the location of cracking. No cracking 
was evident in Figures 3-6 and 3-7. It would have been extremely beneficial if cracking could 


have been observed in Figure 3-6, because the stick and slip zones are readily identifiable. As 
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Table 3.9: Comparison of the analytical versus actual values for the fretting geometry, Al 
7075-T6 (* indicates value could not be determined). 
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Figure 3-7: Picture of the left side of the specimen used for test 33 with the approximate size 
of the stick zone shown (Al 7075-T6). 


Table 3.9 shows, cracking was found on eight of the fretting pads. The cracking found on the 
left fretting pad for test 32 is shown in Figure 3-8. The cracking is on the right side of the 
picture and appears to have started approximately at the edge of contact. The stick and slip 
zones cannot be distinguished, so the exact location of crack initiation cannot be determined 
from this figure. It is interesting that cracking was found on the fretting pads used for test 
32. but no evidence of cracking can be found on the pads or specimen used for test 33. The 
two tests had exactly the same contact loading and the difference in the apphed bulk stress was 
only 7 MPa. It would seem logical that fretting would have initiated cracking in both tests 
and an arrested crack would be produced in test 33 due to the lower bulk stress, but this was 


not found. 


3.4.2  Ti-GAl-AV 


The results of the observations performed on the Ti-6Al-4V test specimens are shown in Table 
3.10. Again it was very difhcult to observe the fretting geometry for the specimens which 


were cycled to failure: however, test 5 provided an interesting specimen. Even though the 
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l'igure 3-8: Picture of cracking on the left fretting pad used for test 32 with the approximate 
circular area of contact shown (Al 7075- T6). 


specimen was cycled to failure, the fracture surface went above the scar on the left side of the 
specimen and preserved the scar. The ratio of c/a was calculated using a coefficient of friction, 
jt = 0.5, which produced fairly good agreement between the analytical and calculated values 
of c/a. The predicted values for the radius of contact were also in reasonably good agreement 
with the actual values observed for the tests. 

The scanning laser microscope was again used to produce photos of the pads and specimens. 
The left fretting pad used in test 5 is shown in Figure 3-9. The shp and stick zones are clearly 
visible in this figure, which is unusual for a test that went to specimen failure. The left side 
of the specimen used in test 5 is shown in Figure 3-10. This is the specimen discussed earlier 
where the fracture plane was above the fretting scar on the left side of the specimen. 

It is interesting that cracking cannot be observed in Figure 3-10. It seems reasonable 
that cracking would initiate on both sides of the specimen since both sides are under the same 
loading conditions. Cracking was observed in the fretting pads of test 2. The cracking in the 
right fretting pad used in test 2 is shown in Figure 3-11. Cracking can be seen on the left edge 


of the scar. Cracking was also found on the opposite side of the scar. but it cannot be seen in 
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Figure 3-9: Picture of the left fretting pad used for test 5 (Ti-GAI-4V). 
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Figure 3-10: Picture of the left side of the specimen used for test 5 (Ti-6Al-4V). 


this photo. The crack appears to have begun near the edge of contact; however, the slip and 


stick zones are not distinguishable and the exact location of crack Initiation is not certain. 
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Figure 3-11: Picture of cracking in the right fretting pad used for test 2 (Ti-6AI1-4V). 
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Chapter 4 


Fretting Fatigue Life Prediction 


Analysis 


4.1 Introduction 


There are two basic approaches used to predict component hfe for plain fatigue: the damage 
tolerant approach and the total life approach. The damage tolerant approach assumes an initial 
flaw exists in the material and the life of the structure is the number of loading cycles it takes 
to propagate an initial flaw to a critical length where failure occurs. In this method only crack 
propagation ts considered. The total hfe approach considers the whole failure crack formation 
process, from initiation to component failure. This work uses both approaches to analyze the 
fretting fatigue phenomenon to provide as much insight into the problem as possible. 

Developing a life prediction method to be used in fretting fatigue apphcations is difficult 
because the effects of fretting are most damaging during the stages of fatigue crack growth 
which have the least developed methodologies: the crack initiation stage and the short crack 
propagation stage. Methods do not even exist for plain fatigue to determine the number of 
cvcles required to initiate a crack and propagate 1t the distance required to use long crack 
propagation techniques. The varving stress fields introduced by fretting make this task even 
more difficult. 

The deleterious effects of fretting on fatigue were shown in the preceding chapter. The 


apphed bulk stress needed to cause fatigue failure in materials under fretting conditions can 
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Figure 4-1: Comparison of the S-N curves for plain fatigue and fretting fatigue for Al 7075-T6 
0 331 MBA): 


be much lower than the endurance strength of the material. The endurance strength of Al 
7075-T6 is 176 MPa; however, under fretting conditions failures were observed at bulk stresses 
as low as 74 MPa. The damaging effects of fretting on fatigue life is shown graphically in 
Figure 4-1 for Al 7075- T6. The constant contact loading series of tests (P = 20 N, Q —15 N) 
has been used for the fretting fatigue S-N curve shown in the figure. Therefore, the fretting 
S-N curve represents only one contact loading condition. The curve could change significantly 
for a different set of contact load parameters. Tlus will be discussed in more detail later in 
this chapter. 

The fretting S-N curve has been shifted down considerably relative to the plain fatigue S-N 
curve in Figure 4-1. Therefore, cracks are initiating and growing to failure in fretting conditions 
where there was no failure for the plain fatigue conditions, even at higher alternating stresses, 
which demonstrates the critical role of fretting in crack initiation. The fretting fatigue curve 
also exhibits an endurance limit. The existence of an “endurance” type limit for fretting fatigue 


is discussed in more detail in the following section. 
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4.2 Total Life Approach 


The total life approach is further subdivided into two categories: tlıe stress-life method and 
the strain-life method. The stress-life method is used in high cycle fatigue (HCF) applications 
where the applied stresses are maintained m the elastic range. ‘The strain-life method is used 
for low cvcle fatigue (LCF) applications where plastic strains make up a significant portion of 
the total strains in the material. The fretting fatigue phenomenon is a HCF problem and the 
stresses calculated in the previous chapter were found to be within the elastic region, so the 
stress-life method has been chosen for further investigation. The stress-life approach is still 
widely used in industry today. 

The stress-life approach is based upon the work of Wöhler performed in the 1850's and 
1860's, who developed the S — N diagram. The S — N diagram is a plot of the alternating 
stress, S, versus the number of cycles to failure, N. The work of Wöhler clearly showed 
a relationship between the applied stress amplitude and the number of cycles to failure and 
indicated an existence of an endurance limit for some materials. The endurance limit is the 
applied alternating stress below which fatigue failure will not occur. Some materials do not 
exhibit an endurance limit. For these materials the endurance strength has been defined as 
the stress amplitude which produces specimen failure at 5 x 10° cycles, or more often shown as 
2N = 1 x 10° reversals. 

A mathematical expression for the stress-life relationship found by Wohler was not produced 
until 1910 by Basquin [26]. Basquin noted that a linear relation existed between the applied 
stress amplitude and the number of reversals to failure when plotted on a log-log scale. The 
relationship is given bv: 


where c, 1s tlie applied stress amplitude, o is the fatigue strength coefficient, bis the fatigue 
strength exponent and 2; is the number of reversals to failure. The fatigue strength coeffh- 
cient, GE is approximately equal to the true fracture strength of the material, oer. The fatigue 
strength exponent, b, is generally between -0.05 to -0.12. This equation has been used quite 
extensively for many years. ‘The values for 6) and b are known for manv different materials 


and can be found in many text books, including references [24] and [22]. 





The Basquin relationship given in Equation 4.1 can be rearranged to give 
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A modification to the equation above is proposed for fretted components such that the number 


of cvcles to failure is now given bv: 
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(1.3) 

where F is the fretting correction factor, which is used to account for the premature initiation of 

cracking in fretted components. ‘To be consistent with earlier notation, ao, has been substituted 

in for og. Equation 4.3 1s very simple and if its validity can be proved, it would be a valuable 

design tool. The critical parameter in the modified Basquin expression is the fretting correction 

factor. F. The values for the fatigue strength exponent. bl and the fatigue strength coefficient. 
/ 


c , are taken to be the same as for plain fatigue. They have not been modified to account for 


fretting. Equation 4.3 can be rearranged such that F is given by: 
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The experimental results given in Table 3.1 have been used to determine /° for each test. The 
analysis of IF is separated into two categories; the constant applied bulk stress results and the 
constant contact stress results. This was done to determine the effect of the contact and bulk 
stresses on IF. F is assumed to be a function only of the bulk and contact stresses. The 
constant bulk stress results are given in Table 4.1. The values of F for the unfailed specimers 
were determined using the number of cycles where the tests were stopped: therefore, these 
values of F are maximum values for the given contact conditions. 

Figure 4-2 1s a graphical representation of the information displaved in Table 4.1. This 
hgure shows the general trend of the data, as the effective stress due to contact is increased. 
the fretting correction factor, F. also increases. The plot includes the unfailed specimens in 


order to display the behavior of F at low contact stresses. The F values for the unfailed 
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449, 500 337 ) 
361.000 366 
530.000 309 
3,300.000 (No failure) 264 
10,680,000 (No failure) 248 
803.000 292 
2,940,000 266 
3,450,000 (No failure) 264 
616,500 285 
2,608.000 (No failure) 280 





049.000 337 


049,000 
3.050.000 


390,000 
258,063 
181,118 

2,197,226 
449,500 





Table 4.2: Fretting correction factors for different bulk stresses, constant contact loads. 


specimens are only approximations, however they do represent the maximum value that the 
fretting correction factor could be for these tests. As can be seen in Figure 4-2, the initial 
slope of the F curve is fairly flat then rises sharply before trailing off again at the lugl contact 
stress values. The initial slope may even be less, because the E values for the unfailed specimens 
are their maximum values. ‘These values may actually be much lower. It must be stressed that 
these results are for one value of op», a different bulk loading will yield a different curve. 

The values of F for the constant contact loading tests are given in Table 4-3. Only the 
tests that were cvcled to failure were included in this table. Table 4.2 has also included test 
number one. The contact loading for this test is nearly identical to the others. ? = 20.8 N 
versus 20.0 N used in the other tests. The plot of IF versus o» is displayed in [igure 4-3 . 

From Figure 4-3, it can be seen that for a given contact stress, /" decreases as the bulk 
loading is increased except for the two low bulk load cases (tests 28 &32). These two cases 


will be discussed again later. 
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Figure 4-2: Plot of the Fretting Correction Factor, F versus the eflective stress due to contact, 
o = 83 MPa. 
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Figure 4-3: Plot of F versus applied bulk stress, while maintaing constant contact loads for Al 
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Figure 4-4: Plot of the combination of bulk and contact loading above which fretting fatigue 
failure will occur for Al 7075-T6. 


The constant contact stress analysis and the constant bulk stress analysis will now be 
combined to provide an estimate of the general behavior of F for various contact and bulk 
loadings. The knock-down factor idea discussed earlier can aid in this effort. ‘Two thresholds 
have been found for the two series of tests. Given o, = 83 MPa, specimen failure will not occur 
il Gee < 201°MIPa, and given oe, = 337 NIPa, specimen failure will not occur if u = (0 MPa, 
These two thresholds can be combined with the endurance strength of the material to produce 
the failure envelope curve shown in Figure 4-4. Specimen failure is predicted for combinations 
of contact and bulk loadings which fall above the curve. 

Figure 4-4 can now be used to analyze Figure 4-2. For occ < 264 MPa fretting will not 
cause specimen failure if o, < 83 MPa; therefore F = 1.0 for these loading conditions. Figure 
4-2 can now be redrawn as seen in Figure 4-5. 

The behavior of the curve in Figure 4-5 is very Interesting. Fretting has no effect on the 
fatigue life of the specimen for a large range of a... However, once the threshold value of Cec 
is crossed, the impact on the fretting correction factor, F., is quite significant. Below the 
threshold, the contact stresses are not high enough to initiate and grow the crack to the length 


required for long crack propagation by the bulk stress. Once the threshold contact stress is 
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Figure 4-5: Revised plot of F versus dec for Al 7075-T6, op = 83 MPa. 


reached cracks will grow enough to be controlled by the bulk stress, but the rate at which this 
occurs is not greatly effected by how far above the threshold the contact stresses are. The 
steep portion of the I curve appears to be a very narrow region surrounding the threshold. 
The slope of the upper portion of the curve is relatively flat. 

The constant contact loading data in Table 4.2 and Figure 4-4 can now be used to estimate 
the behavior of / for other bulk loading conditions as shown in Figure 4-6. The threshold 
values for op = 100 MPa and o, = 125 MPa were determined using the failure envelope curve. 
The shape of the curves were assumed to be similar to that found for o, = 83 MPa. It can be 
seen that the /" curves for the 100 and 125 MPa bulk loading cases pass through the previously 
determined values of 2.38 and 2.0 at Cee = 337 MPa. 

Tests 28 & 32 discussed earlier provide results tliat support tlie validity of tlie curves pre- 
sented above. The F value found for these two cases fell below the curve for the results plotted 
in Figure +3. The threshold Cece determined from Figure 4-4 is 311 MPa for op = 74 MPa. 
Therefore, the contact stress of 337 MPa is approaching the threshold value and could be on 


the steep portion of the F curve, which would result in an F value lower than tle expected 
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Figure 4-6: Plot of the predicted behavior of Jh versus gece for various bulk loadings (Al 7075-T6). 


peak value. 

For a given set of loading conditions Figure 4-6 can be used to determine F, which can then 
be substituted back into equation 4.3 to determine the component life. However, Figure 4-6 is 
only a rough approximation of the behavior of F for various loading conditions. Many more 
tests need to be conducted to determine a better estimate for the behavior of F. The threshold 
contact stress values are the most critical values because it is there where the correction factor 
is changing significantly and hence where fretting fatigue life is changing the most. 

The advantages of this method are that it is simple in form and utilizes the familiar stress- 
life approach widely used in industry today. However, to implement this method will require 
considerable experimental work for many different materials. The method also does not provide 
information on the location of crack initiation. The method uses the contact stresses at the 
edge of contact, but this is not alwavs the location of crack initiation. 


The analvsis performed in this section has produced some very interesting results which are: 


e A rough approximation for the component failure envelope as a function of the bulk and 


contact stresses has been determined and is shown in Figure 4-4. 
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Figure 4-7: Plot of F versus ce. for Ti-GAI-4V (o5 7 300 MPa). 


e The slope of the top of the F curve 1s relativelv flat indicating that the fatigue life 1s not 


greatlv affected bv changes in tlie effective stress due to contact in this region. 


e Ll acts almost as a step function making the determination of the threshold stress values 


of great importance. 


The same analysis was also performed for the Ti-6AI-1V tests recorded in the previous 
chapter. Tests 5-9 from Table 3.7 were used to determine J“ as a function of Cee for op = 300 
MPa. The results are shown in Figure 4-7. This curve 1s based on limited experimental data, 
but the same type of behavior as seen with the aluminum alloy seems to be present for the 
titanium alloy. For the given bulk loading, / = 1 until ges Z 390 MPa then it quicklv jumps 


fore er. 


4.3 Damage Tolerant Approach 


4.3.1 Al 7075-6 


The second life prediction approach used to analyze the fretting fatigue problem is the damage 


tolerant approach. Using the experimental results from the preceding chapter. this work 


Al 


attempts to determine the relative magnitudes of the crack nucleation lifetime aid the long crack 
propagation lifetime. The crack nucleation lifetime used in this work is the hfe spent initiating 
and growing the fatigue crack to the length where the fretting stresses are approximately zero. 
This approach is used by Szokwinski and Farris [30] in their fretting fatigue life prediction 
method. The long crack propagation hfe then is only influenced by the bulk stress. Linear 
elastic fracture mechanics can be used to determine this life. This long crack propagation 
life can then be compared with the total fatigue hfe from the experiments to determine the 
relative magnitudes of the two stages of fatigue life. The approach is demonstrated using the 
A] 7075- T6 experiments and is applied in the following section to the T1-6A]-4V experiments. 

Paris, Gomez and Anderson [31] in 1961 and Paris and Erdogan[32] in 1963 demonstrated 
that a power law relationship exists between the crack growth increment da/dN and the stress 


intensity factor range AA. This relationship is given by: 
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where C' and m are scaling constants, and 
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Observations were made of the failed test specimens to determine the shape of the propagating 
crack. It was determined that the cracks could be best characterized as surface elliptical cracks 
with a ratio of a/c = 0.9. A diagram of this geometry is shown in Figure 4-8, with the thickness 
and width of the fatigue specimens also given. 

The stress intensity factor equation for the surface elliptical crack geometry was determined 
by Newman and Raju [33] in 1981. For a finite plate subjected to a remote uniform tensile 


stress, S, the stress intensitv factor is given by: 
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Figure 4-8: Diagram showing the elliptical crack geometry. 


where 
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The expression for F was obtained using a systematic curve-fitting procedure [33]. Note this 
F is not the same F described in the previous section and a in these expressions refers to the 


crack geometry and not the contact radius. The expression for F is given by: 
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Table 4.3: Values used in the Paris life calculations for Al 7075-T6. Consistent units for C are 
meters and MPa. 
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The function, fw is the finite width correction. 

The above equations can now be used to find A given any a. The expression for the stress 
intensity factor is now substituted back into the Paris form, Equation 4.5 and the equation Is 
numerically integrated from an initial crack length ag to the critical crack length ay to find the 


crack propagation life. Using @ = 90°, the propagation life, Ny» is given by: 
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C(Ae) 'Jw (9.67) a? [1.05 + 0.27 (2)? - 0.145 (2)%| [sec (102 /2)] 


The values used in the propagation life calculations are displayed in Table 4.3 [35]. The 
value used for ag is approximately four times the contact radius for the tests conducted on the 
aluminum alloy. The stresses due to contact decay away very quickly and are approximately 
zero at a depth of 1.5a. However, if 1.5a is used for ag, the calculated propagation life is greater 
than the experimental life, which cannot be true. ‘Therefore, ag was set equal to da. The 
value used for ag must also be long enough. such that it includes at least a few grain diameters. 
The average grain size for the Al 7075- T6 specimens was 60 jum [23], so the ag used in these 
calculations encompasses several grain diameters. As can be seen from Table 4.3, only one-half 
the thickness of the specimen was used in the calculations. This is due to the presence of two 
cracks in the specimens. When the fracture surfaces of the failed specimens were analyzed, 
two cracks were observed. One on each side of the fatigue specimen. 


Using Equation 4.14. the long crack propagation life was calculated for each test, using o, 


AA 


for Ac. The stress intensity factor range was also calculated for each test using Equation 4.6, 
at a = ag, and S = Cp to examine if AK > Ag, the threshold stress intensity factor range 
at this crack length. The AKo for Al 7075-T6 at & = —1 appears to be approximately 2-3 
MPaym. This value was determined from a plot in [34]. The plot did not contain R = —1 
crack growth data: however, using the displayed trends for the other It ratios an estimate for 
the AK for R = —1 was made. The calculated values of the propagation life. Ny, and the 
stress intensity factor range at a = ag for each test are given in Table 4.41. 

As can be seen from Table 4.4, the long crack propagation life is a verv significant portion 
of the fretting fatigue life for many of the failed specimens. For these specimens, a verv small 
portion of the life was used to nucleate cracks. The calculated AK values in Table 4.4 provide 
insight into some of the tests where component failure did not occur. It can be seen that for 
Tests 14-17 and Tests 21-24, AK « 2.5 MPaym for each of these tests, which is very close to 


the threshold value for this material. 


4.3.2 Ti-6Al-4V 


The method used for Al 7075-T6 above is applied to Ti-6A]-1V in this section. The crack 
geometry for the titanium alloy specimens was very similar to that found above. Again elliptical 
surface cracks were found; however, cracking was found on only one side of the specimen for this 
material. The values used to calculate the propagation life for this material are shown in Table 
4.5 (a/c corresponds to the elliptic crack geometry). The values used for ag are approximately 
1.5a and da. The first value of ag — 1.5a, corresponds to the depth at which the tensile stresses 
due to contact are approximately zero. The second value of ao — 4a is used to be consistent 
with the calculations performed earlier for A1 7075- T6. 

The values for m and C were determined using curves found in [37] and [36]. The fatigue 
crack growth curve in [36] was for solution treated and over-aged (STOA), Bi-modal Ti-6Al-4V. 
while that found in [37] was for Bi-modal Ti-6Al-4V with no heat treatment. Both curves were 
made using R = 0.1. The values found for C and m for the two cases are shown in Table 4.6. 
The experimental results displayed in Chapter 3 are for Bi-modal Ti-6Al-4V with no special 
heat treatment, at R= —1. The value of m for this section was set to 5.0, since the applied 


stress ratio has a small effect on m. The value of C was determined by approximating the 









































































Experimental Ny (Cycles) | Nj, (Cycles) | AK PN 




































449,000 408,000 
480,000 371.000 
395,000 371,000 
361.000 408,000 
551.000 991,000 
530.000 408,000 
3,300,000* 108,000 
10,680,000* 103,000 
803,000 408,000 
2.940.000 408,000 
3,450,000* 108,000 
616,500 408,000 
1,410,000" 855,000 
2,186,000" 2,280,000 
2.500,000* 1.400,000 
16 2.666,000* 1,400,000 
I 2 005.000" 1,400,000 
18 2.608 ,000* 408,000 
19 949.000 408,000 
20 516,000 807,000 
21 1.540,000* 1.970,000 
22 2.940.000* 1.230.000 
23 1.77 7.000* 1.600.000 
24 297,000 3.000 
25 135,000 3,650 
26 480,000 6,920 
2 919,000 13,400 
28 3.050.000 612.000 
29 390.000 194,000 
30 258.063 19,300 
3 1211 38,300 
32 2,497,226 646,000 
35 3,280,000* 961,500 2^ 





Table 4.4: Long crack propagation life for Al 7075-T6. 
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1.0 
5.0 mm 
5.0 mm 
5.0) 
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200 um, 180 uum 
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Table 4.5: Values used for the propagation life calculations for Ti-6AI-4V. 


| | STOA, Bi-modal Ti-6AI-4V | Bi-modal Ti-GAI-4V 
C 8.13x 107 Ps 
m 9:12 5.0 


Table 4.6: Values of C and m for STOA, Bi-modal Ti-6Al-4V and Bi-modal Ti-6AF4V (units 
are meters and MPa). 






effect of R on C. 

Using the values displayed in Table 4.5, the propagation life was calculated for the nine 
tests shown in Table 3.7. "The propagation lives corresponding to the two different values of 
ag are shown in Table 4.7. As can be seen from the table. the propagation hfe correlates very 
wel with the actual experimental life of the failed specimens when ag = 1.5a is used. Using 
ag = 1.5a, the difference in fatigue Hfe for the various tests can be attributed almost entirelv to 
the difference in propagation life. Using ag = 4a results in much lower values of Ny,. ‘Table 
3.7 also displavs the calculated AK values at ag = 1.5a. The threshold stress intensity factor 


range, AKo, is approximately 3-7 MPaym [37]. 


RED AK (MPa itas c 150) 
691,612 1,010,000 134,000 de 
109.514 06.800 1,040 
107,879 61,200 8,130 
31.540 (interrupted) 100,000 13.300 


304.456 247,000 32,900 

Kee 239,000 31.800 

143,683 156,000 20,700 6.96 

258.916 239,000 31,800 Be 
3,257.000(no failure) 239,000 31.800 6.38 





Table 4.7: Ti-6Al-4V long crack propagation life. 





Chapter 5 


Palliatives 


5.1 Introduction 


The work performed in the previous chapters provide the necessary background for a discussion 
of fretting fatigue palliatives. The methods developed earlier provide the criteria that can be 
used to determine how and whv certain palliatives work while others do not. 

This work will focus on two fretting fatigue palliatives: shot peening and coatings. Both 
have been shown to be beneficial in increasing fretting fatigue life. Chivers and Gordelier [41] 


provide a detailed study of several possible palliatives including shot peening and coatings. 


5.2 Shot Peening 


Shot peening has been used for many years as a way to improve the fatigue properties of 
components in plain fatigue. The benefits of shot peening were found quite bv accident bv 
the Buick Motor Division of General Motors Corporation in the late 1920's [42]. It was found 
that valve springs that had been blasted with grit to remove scale had greatly improved fatigue 
properties. It was first thought that the benefits derived from shot peening were due to the 
work hardening of the surface. However, it was soon discovered that the benefits were caused 
bv the residual stresses introduced at the surface due to shot peening. 

In the shot peening process, relatively hard particles are shot against the metal being treated. 


The particles are generally spherical in shape and are given enough velocity such that the surface 





= Fretting & Plain Fatigue (shot peened) 


Plain Fatigue (not peened) 
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Fretting Fatigue (not peened) 





Cycles to Failure 


Figure 5-1: Plot showing approximate relationship between shot peened and unpeened speci- 
mens for plain and fretting fatigue. 


of the metal is indented during impact. The indentations occur due to local plastic vielding and 
as the surface material attempts to expand, it is restrained by surrounding material below the 
surface which was not deformed plastically. This restraint of the plastically deformed surface 
laver causes it to be compressively strained and is residually stressed in compression. 

There have been numerous studies conducted showing the benefits of shot peening for 
plain fatigue conditions. Shot peening has also been demonstrated to be effective in fretting 
fatigue conditions as well, which is discussed in [38], [39] and [41]. In many materials, shot 
peening lias been shown to increase the fatigue strength such that the fretting fatigue strength 
iS approximately tlie same as the plain fatigue strength. <A simple plot displaying the typical 
relationship between shot peened and unpeened specimens for fretting and plain fatigue is 
shown in Figure 5-1. This figure is based upon the work of Bignonnet [39], which studied the 
effects of shot peening on fretting fatigue. However, this figure is not true for all materials. It 
is only meant to display the benefits typically observed in shot peened materials. 

As discussed earlier, the benefits of shot peening are derived from the residual compressive 
stresses introduced into the material during the process. ‘The maximum value of the residual 
compressive stress, often called the magnitude, is not affected by the shot peening process as 
long as the shot used is at least as hard or harder than the material being treated [43]. The 
magnitude of the residual compressive stress 1s material dependent, but has a value of at least 


one-half the vield strength for that material [43]. The depth of the compressive layer is the 





other important parameter in shot peening. The depth of the compressive layer is a function 
of the hardness and ductility of tlie peened material, the state of strain in the material when 
the peening was performed and the characteristics of the shot and shot stream [42]. A typical 
residual stress profile due to shot peening and the residual stress profile used in tliis analysis are 
shown in Figure 5-2. The residual stress profile used in this analysis assumes a flat maximum 
compressive stress, which goes directly to zero at a given depth. ‘The assumed residual stress 
profile should actually have a tensile region to balance the residual compressive stresses to 
satisfy equilibrium. This is shown in Figure 5-3, where 07. is the residual compressive stress 
and or is the residual tensile stress required for equilibrium. For equilibrium to be satisfied 


area Al must equal area 42. Therefore the area balance is given by: 


(Odi ae) =e) (ary |. (5.1) 


and the ratio of the residual stresses is given by: 


i y a 

C 
If (t — pd) > (pd) then 07; — 0. It is assumed in this work that (t — pd) > (pd), therefore the 
residual tensile stress has been neglected. Note that the actual maximum compressive stress 
occurs at some depth below the surface. The actual depth of residual compressive stresses for 
various materials is given in [43] as well as some residual stress profiles. 

The benefits derived from shot peening can also become detrimental to the material if the 
residual compressive stresses are somehow relieved. The process which caused the residual 
stresses to be induced into the material also damages the surface of the material and may even 
introduce microcracks. These microcracks could drastically reduce fatigue hfe if the residual 
compressive stresses are reHeved. This has been shown to be the case in both plain fatigue 
and fretting fatigue. Gray, Wagner and Liitjering [44] show the detrimental effects on plain 
fatigue life of shot peening when the residual compressive stresses are thermally relieved. The 
fatigue strength of the material is significantly reduced for the shot peened, stress relieved 
specimens as compared to the original electrolytically polished specimens. Figure 5-4. from 


[44]. demonstrates this effect. In Figure 5-4, SP is shot peened, EP is electrolytically pohshed 
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Figure 5-2: Typical stress profile after shot peening and the assumed stress profile used in this 


work. 


pd = peening depth 
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Figure 5-3: Residual stress equilibrium diagram. 
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Figure 5-4: Diagram showing the detrimental effects of shot peening if the residual stresses are 
relieved [44]. 


and SR is stress relieved. [rom this figure it can be seen that the fatigue strength for the 
shot peened and stress relieved specimens is about one-half that of the electrolytically polished 
specimens at 1x10^ cvcles. the reduction in fatigue strength is quite significant and certainly 
cannot be ignored. The shot peened and stress relieved curve shifts down relative to the 
electrolvticallv polished specimen curve indicating that the difference between the two curves 
is in the crack initiation stage, just as was found for fretting. The material used in [44] was 
Ti-6A1-4V with a fine lamellar microstructure and the tenıperature was 20 °C for this set of 
test results. The residual stresses were relieved thermally using a heat treatment of 600 °C for 
one hour. These results were for plain fatigue experiments. Bignonnet [39] displays similar 
results for fretted specimens. The fretting fatigue strength of the shot peened, stress relieved 
specimens was similar to that of the unpeened specimens. Waterhouse [38] found similar results 
in his experimental work. These works demonstrate the importance of the compressive residual 
stresses and the importance of not allowing them to be relieved. 

Laser shot peening does not have the shortcomings discussed above for the typical shot 


peening process. In the laser shot peening process the surface is unchanged by the peening 


A) 





process. Shock waves are introduced into the material that are capable of plastically vielding 
the material and producing compressive residual stresses. 

The residual stresses can be relieved either mechanically or thermally. The approximate 
temperature at which the residual compressive stresses begin to dissipate for various materials 
can be found in [43]. The residual stresses are relieved mechanically if the total stresses in the 
material exceed the dynamic elastic limit of the material. ‘The total compressive stress is made 
up of the shot peening residual stress and the stress due to external loading. The residual 
stress is relieved by the amount that the total compressive stress exceeds the dynamic elastic 
limit of the material. 

The objective of this section is find the optimum peening depth for various fretting fa- 
tigue loading conditions. The magnitude of the residual compressive stresses are fixed by the 
material, therefore the designer has no control over this aspect of the shot peening process. 
However. the depth of residual stress can be influenced by variations in peening parameters as 
shown in [43]. Therefore, the depth of peening becomes the variable design parameter that 


can be optimized by the designers. 


5.2.1 Material Yielding Analysis 


A systematic analysis of the shot peening process applied to the sphere on flat surface contact 
geometry has been conducted in this work. The material used for this analysis is Ti-6Al- 
4V. The pertinent material properties for Ti-6Al-4V are listed in Table 3.6. The first task 
was to find the conditions when the residual compressive stresses would begin to be relieved. 
According to [43]. the compressive stresses begin to dissipate thermally at 480 C for Ti-6AI-4V. 
Mechanical stress relieving begins to occur when the material stresses exceeds the elastic Hmit 
of the material. The stress fields for a sphere fretting against a flat surface have already been 
discussed in Chapter 2 of this work. The surface stresses for this geometrv for Al7075-T6 are 
displayed in Figure 2-2. "The stress profile for Ti-6A]-4V is almost identical, the magnitudes 
are different, but the shape of the curves are the same. 

The analvsis performed in this chapter assumes the onset of global sliding. This condition 


provides the highest stresses possible for a given contact loading. If global sliding is assumed, 


the normal load. P, becomes the critical contact parameter since the tangential load, Q, will 
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always be equal to uP. The loading which causes the total effective stress, Octor, to become 
equal to the yield strength is when material yielding is assumed to occur and is sliown in Figure 
5-5. The data used to make this figure were determined using a MATLAB script file. <A 
MATLAB routine was written which first scans the surface of the material, then incrementally 
scans below the contact surface calculating the stresses at each increment. producing a matrix 
of stresses. The surface is scanned along the y = 0 hne. where maximuni material stresses are. 
The increment used in the x-direction is 0.02a. The increment in the z-direction is 0.025a. 
The program then finds the maximum effective stress from the matrix of stresses determined 
during the scan. The stresses at the contact edges are calculated separately. The maximum 
stress is then compared to the vield stress. If the maximum effective stress is less than the 
vield stress, the contact loading is increased and the stress calculations are performed again. 
The process is continued until the contact loading which produces yielding 1s determined. This 
was done for five different coefficients of friction and eight different bulk stress conditions. The 
stresses at the edge of contact are calculated separately because the stresses are indeterminate 
at the contact origin (r — y — z — 0) and the origin would be included in the surface scan if 
both contact edges are also included. The program used for this analvsis is shown in Appendix 
A. 

In Figure 5-5 the vielding curves have been plotted. which show for a given bulk stress, 
cy ànd a given coefficient of friction, u, the minimum contact pressure. Fo, that is needed to 
cause the total effective stress to exceed the yield stress of the material. For example if u = 0.7 
and o, = 0.loy then the maximum contact pressure which produces yielding is ~0.280y. This 
example is shown in Figure 5-5. The maximum contact pressure is normalized bv the yield 
strength of the material in this plot. A residual stress of -0.850y is used for this study. This 
value for the residual compressive stress was determined using an actual shot peening stress 
profile from [44] and a graph of maximum compressive residual stress versus tensile strength in 
[13]. 

As can be seen in Figure 5-5, the contact loading which causes material yielding increases 
as the applied bulk stress increases until 0, = 0.doy. The leading edge of contact, where the 
compressive stresses due to contact are very high, is the location of vielding for applied bulk 


stresses less than one-half the vield strength. The location of yielding shifts to the trailing 
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Figure 5-5: The combination of contact and bulk loading which causes the total effective stress 
to equal the vield stress plotted versus the coeflicient of friction. 


edge, where the tensile stresses due to contact are very high, when op > O.d0y. Therefore, 
vielding is controlled by the compressive stresses until the applied bulk stress reaches a value 
of one-half the vield strength. This explains why the maximum contact pressure which causes 
vielding actually increases with increasing bulk stress when op < 0.1oy. Figure 5-5 also shows 
the effect of friction. As the coefficient of friction increases, the contact pressure needed to 
cause vielding decreases. [t can also be seen that the contact loads needed to cause yielding are 
relatively low for the coefficient of friction typical for Ti-6Al-4V (0.4-0.5); therefore, residual 
stress reheving due to vielding cannot be ignored. Yielding first occurs at the surface for u > 0.3. 


Since vielding usually occurs first at the surface, the depth of vielding must also be found. 


5.2.2 Depth of Yielding 


The depth of yielding must be found to determine to what depth the residual stresses are being 
relieved. The von Mises effective stress is again used to determine the yield conditions, 1.e. 
when Oe > oy the material is considered to have vielded. A MATLAB routine was also 


written to perform this task. The routine used to perform this task is similar to the previous 
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Figure 5-6: Plot of yielding depth. z versus maximum contact pressure for various bulk loadings 
E10). 


one. Again the surface 1s scanned using an increment of 0.02a in the x-direction with y — 0 
for all points. The stresses are then determined below the surface using an increment of 0.02a 
in the z-direction, which results in a matrix of stresses. The stresses at the edge of the contact 
are not calculated separatelv in this routine because the stress profile is much smoother at the 
edges below the contact surface as shown in Figure 2-3. The depth of maximum stress is found 
and this depth became the starting point for the iteration to find the depth of yielding. The 
depth of maximum stress is always the surface for u > 0.3. ‘The scanning process begins again 
at the depth just determined using a reduced increment in the z-direction of 0.0la. The depth 
is increased until the effective stress is no longer greater than the vield strength of the material. 
Figure 5-6 displays the results of this analysis for z = 1.0. The depth. z is normalized by the 
radius of contact, a, the maximum contact pressure, Po. 1s normalized by the vield strength of 
the material, cy. Again the residual stress used in this analvsis Was Oresidua! = —0.850y. The 
MATLAB script file used in this section 1s shown in Appendix B. 

The results displayed in Figure 5-6 follow the expected trends, the depth of vielding increases 


as the contact load is increased and decreases as the bulk loading is increased until 0, = 0.50 y. 
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Figure 5-7: Depth of yielding versus maximum contact pressure for various bulk stresses (ji = 
0.8). 


The depth of vielding decreases with increasing bulk stress due to the large residual compressive 


tress assumed in the analysis. This same analysis was repeated for three more coefficients of 


N 


friction, 44 — 0.8, u = 0.5 and u = 0.3. The results for the three other coefficients of friction 
are shown in Figures 5-7, 5-8 and 5-9 respectively. As can be seen in these figures. the depth 
of vielding is reduced significantly as the coefficient of friction is lowered, which again points 
to the possible benefits of coatings in fretting situations. When the coefficient of friction is 
reduced to 0.3, there is very little vielding even at verv high contact loadings for o, > 0.4. 
The depth of vielding curves developed in this section can be very useful for design, specif- 
ically in the material treatment stage of design. Designers could use these plots with the 
known contact and bulk loadings and the coefficient of friction and very quickly determine 
the minimum depth of shot peening to ensure that the compressive residual stresses are not 
mechanically relieved. For example. if jp = 0.8, Po/o, = 0.8, and o, = OU.loy, the depth 
of vielding, z/a ~ 0.35. Therefore, the minimum depth of shot peening to avoid mechanical 
relieving of the residual stresses is z/a > 0.35. This example is shown in Figure 5-7. The 


depth of vielding curves provide important design information, but thev only provide part of 
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Figure 5-8: Depth of yielding versus maxımum contact pressure for various bulk stress loadıngs 
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Figure 5-9: Depth of vielding versus maximum contact pressure for various bulk stress loadings 


(gt = 0.3) 
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the information needed to determine the optimum shot peening depth. These curves provide 
the minimum shot peening depth to avoid residual stress relieving but they do not provide 
information pertaining to the tensile stresses at the crack tip. An additional analvsis is needed 


to provide this information. which is discussed in the next section. 


5.2.3 Maximum Principal Stress Analysis 


The predicted success of shot peening is based upon two basic assumptions. First, a crack will 
not propagate unless a tensile stress forces it open near tlie tip [43]. The second assumption 
is that a crack tip will not open while a compressive force is acting upon it [43]. These 
assumptions provide the criteria that is used to determine the optimum shot peening depth. 
Using the two assumptions stated above, the depth of peening to prevent crack propagation 
should be the depth at which the maximum principal stress, 0; < 0. Therefore, inside this 
laver of material affected bv shot peening, no region of tensile stresses will exist. This does 
not mean that the entire depth of the material will have no tensile stresses. ‘This only apphes 
in the shot peened region, but that is enough to stop the crack from propagating. 

A MATLAB program similar to the two discussed above was written and is shown in 
Appendix C. The surface is scanned along the x-direction determining the stresses at each 
point using an increment in the x-direction of 0.02a, again y = 0 for all points. At each point 
the maximum principal stress, 71, 1s determined creating a vector of maximum principal stresses. 
The maximum of this vector is then found and compared to zero. If the maximum of the vector 
is greater than zero. then the depth is increased by 0.01a and the process continued until the 
maximum principal is less than or equal to zero. This depth then became the required shot 
peening depth. The shot peening depth, z, which produced a maximum principal stress less 
than or equal to zero was found for many different loading conditions. Just as in the vielding 
depth curves, the depth is plotted versus maximum contact pressure for several different values 
of applied bulk stress. Figure 5-10 displavs the depth z. at which 01 < 0, for 4 =1.0. Again 
the depth has been normalized by the contact radius, a and the maximum contact pressure, 
Po. is normalized by the vield strength of the material. This analysis was repeated for three 
other coefhcients ol friction. 1 = OS 1 كك‎ 0) ance) TNE CES US 


are shown in Figures 5-11 and 5-12 respectively. The curves for u = 0.3 are not displayed 
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Figure 5-10: The depth at which the maximum principle stress is less than or equal to zero 
plotted versus maximum contact pressure. for magnitudes of bulk stress (ji — 1.0). 


because the required peening depth was less than for 0.10a for even the highest load conditions 
analvzed. The plots again show the advantages of lowering the coefficient of friction. There is 
a significant reduction 1n z/a as the coefficient of friction is reduced. As expected the required 
shot peening depth increases as the applied bulk stress, o; is increased. As seen in the figures, 
the increase in peening depth with increasing a, becomes quite significant when cp > V.Joy. 
The plots generated in this section cannot be used in isolation, they must be integrated with 
the plots from the previous two sections of this analvsis. ‘Together, these sections provide the 
designer the information necessary to determine the optimal shot peening depth for the given 


loading and material conditions. 


5.2.4 Integrated Yielding and Maximum Principal Stress Analysis (Opti- 
mum Shot Peening Depth) 


The depths determined from the two preceding sections provide two different threshold val- 
ues. The depths output from the material vielding analysis provide the minimum depth to 


which the specimen must be peened to prevent mechanical relieving of the residual compres- 
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Figure 5-11: The depth at which the maximum principal stress is less than or equal to zero 
plotted versus maximum contact pressure for various magnitudes of bulk stress (u = 0.8). 
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Figure 5-12: The depth at which the maximum principal stress is less than or equal to zero 
plotted versus the maximum contact pressure for various magnitudes of bulk stress (jz = 0.5). 





sive stresses. The depth of vielding analysis determines the minimum peening depth to ensure 
that the residual compressive stresses remain in the material. The maximum principal stress 
analysis assumes the residual compressive stresses are still present and determines the peering 
depth required to prevent crack growth. Therefore, the depth of yielding analysis provides a 
minimum shot peening depth, but does not provide all information necessary to determine the 
optimal peening depth. The remaining information is provided bv the maximum principal 
stress analysis, which determines the peening depths at which a; <0. This analysis assumes 
the crack will not grow in a region of pure compressive stresses. ‘The two analysis must be 
combined to form the conıplete picture, which has been done in Figure 5-13 for u = 1.0 and 
oy = O.4oy. As can be seen from Figure 5-13 the limiting curve depends upon the contact 
loading, a cross-over point occurs at (Po/oy) = 1.15. If (Po/ay) S 1.15 the maximum prin- 
cipal stress analvsis is the limiting case and controls the optimum peening depth, however if 
(Po/oy) 2 1.15 the depth of vielding analvsis becomes limiting and determines the required 
peening depth. This example shows the importance of examining both criteria. Án example 
on how to use the figure is also displayed in Figure 5-13. H (14/0y) =0.8 then the optimum 
peening depth. z/a — 0.4. The designer may choose to shot peen to a slightlv greater depth 
to ensure no tensile stresses exist in the material at that depth. 

A second example is given in Figure 5-14 for u = 0.8 and 0, =0.1l0y. As can be seen in 
this figure. the intersection between the two curves occurs at a much lower contact pressure 
than in the previous example. Contrary to tlıe previous example. the depth of vielding is the 
limiting curve over most of the contact loading space for the given u and a,. These curves can 
be generated quite easilv from the figures provided in the previous two sections for many other 
combinations of u and gp, however only two examples will be displayed in this work. In Figure 
9-14 the optimum peening deptli. 2/a, was again found for a contact load of (Pof/oy) = 0.8. 
In this example the optimum peening depth was found to be; z/a ~ 0.36. The shot peening 
depth found in this example would be the minimun recommended because material yielding is 
still occurring at this depth. The material is just at tlie yield point, but it is recommended 
that the material be shot peened to a slightly greater depth to ensure no material yielding and 


hence. no mechanical residual stress relieving will occur. 


7 















L س سے‎ p MEME = 
T ZA 
EECH rr MN NC MR E 

- Ta ] 

02 = vm 7 

H É | 4 

[ E : | = “max c, analysis |: 

0:1 r^ pO DIS SHE 
- P d i S— yielding analysis 
t 
0 
0.2 0.4 0.6 0.8 1 1.2 1.4 
P (o, 


Figure 5-13: Optimum peening depth plotted versus maximum contact pressure for u = 1.0 
and 0, = 0.loy. 


0.7 


0.6 =- | 一 二 yielding analysis 
— —— max c, analysis 


E 


0.5 


Ep 








一 
Ka - 
- 
E 
0.4 - 一 -- 
® E 
3 E ] 
0.3 7 = 
E 0 m 8 
E anna ae oe 
DEZA kT i EI Ce MCN ME - 
L J 
I p 
= 7 
0.1 C EECH en S = 
0 
0.2 0.4 0.6 0.8 1 1,2 1.4 
Pic 
0 1 


Figure 5-14: Optimum peening depth plotted versus maximum contact pressure for u = 0.8, 
و12 (اكره‎ 





Parameter Value 





Table 5.1: Parameters and respective values for the threshold stress intensity range analvsis. 


5.3. Stress Intensity Range Threshold Analysis 


"The maximum principal stress analysis is a conservative approach. An alternative to the 
maximum principal analvsis is an analysis based on the threshold stress intensity factor range 
AKo. At stress intensity factor ranges below the threshold. cracks either remain dormant or 
erow at undetectable rates [24]. This could be used in place of the maximum principal stress 
analvsis to determine the optimum peening depth. Using the values listed in Table 5.1, a 
MATLAB script file was written to determine the peening depth such that AA applica « A Ko. 


The expression used for AN applica 15 given by: 


no 
Cl 
CD 

Sy, 


AKappliea = Y (Ao) v Ta 


where Y is a geometrical factor. 

The analysis was performed using op = 0.60y, and u = 0.8. The material stresses in the 
x-direction were used to determine Ao. Figure 5-15 shows the results of the analysis, where 
the required peening depth is plotted versus the maximum contact pressure. In this figure the 
results of the maximum principal stress analvsis for the same loading conditions are also given 
for comparison. As can be seen in the figure, the maximum principal stress analvsis is the 


more conservative analvsis. 


9.4 Protective Coatings 


The protective coatings analysis in this work focuses on the abilitv of the coating to change 
the coefficient of friction of the svstem. Again a sphere on flat plane geometry is used in this 
section. A change in the coefficient of friction has a direct impact on the contact stresses. The 


contact stresses displaved in Equations 2.10 and 2.11 are shown again below in Equations 5.4 
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Figure 5-15: Comparison of the required peening depths for the maximum principal stress 
analvsis and the threshold stress intensity range analysis. 


and 5.5. The stresses at the edge of contact, at the surface, are given by: 








1-2 4 
Cp = 15 = SE San = ame (5.4) 
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Therefore, reducing the coefficient of friction from a value of 1.2 to 1.0 will reduce the stresses 
from a, = 2,160) to oz = 1.822Fo and o, = 0.291 Fo too, = 0,220 F9, assuming Ts 0,5 50 
it can be seen that changing the coefficient of friction has a significant impact on the stresses 
induced by contact and therefore. on the fretting fatigue life of the structure. Figure 5-5 also 
demonstrates the benefit of reducing the coefhcient of friction. It must be remembered that 
the equations shown above were for global sliding conditions. However, Chivers and Gordelier 
in [16] have shown that for yz > 0.2, reducing 4 will also reduce stress for stick-slip conditions. 
No experiments have been conducted as part of this work in which the coefficient of friction 


was purposely altered to examine the actual effect on fretting fatigue life. However. Chivers 


=] 


ZA 





and Gordelier [41] performed a number of experiments for various fretting fatigue palliatives 
including shot peening and protective coatings and their results do support the validity of the 
discussion above. They examined two types of non-metallic coatings. Both were graphite-filled 
resins. one of the resins was a polyimide and the other was an isocyanate epoxy. The improve- 
ments in fretting fatigue behavior produced by the coatings were quite significant and were 
observed at all stress levels. Chivers and Gordelier |-41].define a term called the improvement 


factor, Jy. to measure the success of the various palliatives being investigated. [y is given by: 


1) (5.6) 


UTE 


where gp is the run-out stress under fretting conditions with a palliative, o; is the run-out 
stress under fretting conditions and no palliative and o is the run-out stress without fretting. 
Chivers and Gordelier define the run-out stress as the stress at which failure does not occur 
in about 10° cycles. The improvement factor values were found to be Jy = 0.3 for graphite 
in polyimide resin and Jy = 0.55 for graphite in epoxy resin. For comparison purposes the 
improvement factors found for shot peening were عم[‎ = 0.47 for low intensity peening and 
Iş = 0.43 for high intensity peening. It is interesting that the graphite in epoxy resin coating 
actuallv outperformed shot peening. Chivers and Gordelier also performed experiments using 
a copper shim and the combination of a copper shim and graphite in epoxy. The use of the 
copper shim was intended to reduce the contact stresses at the surface of the material due to 
fretting, which from the previous analysis, again should improve fretting fatigue behavior. ‘The 
improvement factors found were If = 0.63 with copper shim and /f = 0.77 using both a copper 
shim and graphite in epoxy. Flat ended cylindrical pads made of 1Cr-Mo steel were fretted 
against rectangular 3.5Ni-Cr-Mo-V specimens for the experiments conducted by Chivers and 
Gordelier in [41]. The reduction in fatigue life due to fretting is directly related to the stresses 


at tlhe contact surface and if tle stresses can be reduced the fatigue hfe will be improved. 





Chapter 6 


Conclusions and Future Work 


6.1 Conclusions 
The kev findings in this work are summarized below. 


e The surface stresses due to contact can be used to characterize the damaging effects of 
fretting on the fatigue properties of a material. The critical fretting parameters, such 
as: the normal contact load, the oscillating tangential load and the coefficient of friction 
are all incorporated into the expressions for the contact stresses. It was shown that for 
a given bulk loading condition. the fretting fatigue life is reduced as the effective contact 


stresses are increased. 


e A fretting fatigue knock-down factor curve has been produced for Al 7075- T6. It shows 
the approximate relationship between the effective stress due to contact and the magnitude 


of the fretting fatigue knock-down factor, A I’. 


e lt was shown that the use of a single correction factor to modify the Basquin relationship 
for fretting fatigue would require extensive experimental work. The correction factor can 
be expressed as a function of the effective stress due to contact for a given bulk loading 
condition. However, this relationship changes as the bulk loading 1s changed. Therefore, 
families of correction factor curves would need to be developed for various bulk loadings. 


The experimental work required may be tenable if the correction factor was determined 





to be a function of the contact stresses onlv and not dependent on the contact geometry 


for a given material. 


e A failure envelope curve lhas been generated for Al 7075-T6. This curve shows the 


approximate combination of bulk and contact loading that will produce component failure. 


e The damage tolerant approach analvsis has been used to bound the crack initiation and 


short crack propagation stages of fretting fatigue. 


e A new method has been proposed to determine the optimum peening depth for a com- 
ponent that undergoes fretting using the stress analysis provided by Hamilton [14]. The 
method was applied to a sphere on flat surface contact geometry, but is valid for any con- 
tact geometry that has analytic expressions for the stress fields. The optimum peening 
depth has been determined for Ti-6Al-4V for various combinations of bulk and contact 
loadings. for a given coefficient of friction. This analysis was performed for four different 
values of the coefficient of friction, so that the effect of coatings on peening depth could 


be analyzed. 


6.2 Future Work 


More experimental work needs to be performed to determine a more accurate failure envelope 
curve. This could prove to be a very beneficial design tool. especially if it was shown to be 
geometry independent. Theoretical work is also needed for this curve to better understand why 
the given combination of contact and bulk stresses lead to component failure. This will require 
a thorough examination of the crack initiation stage of fretting fatigue crack growth. 

The shot peening analysis provided in this work also needs experimental validation to de- 
termine the actual depth: of mechanical stress relieving due to the contact stresses. If the 
compressive residual stresses due to shot peening are relieved, the damage to the surface by the 
shot peening process may cause a significant reduction in the fatigue strength of the material 
as was discussed earlier. 

More experimental work also needs to be done to determine the exact location of crack 


initiation. The pictures displayed in the previous work do not provide conclusive evidence for 





the location of crack initiation. Interrupted tests are recommended because they could provide 
not only the location of crack initiation, but could also provide information on the number of 
cvcles required for crack initiation. More analysis is also required for the predicted location of 
crack initiation for fretted components. The stress analysis approach does not always predict 


the correct location of crack initiation. 
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Appendix A 


Ti-6Al-4V Yielding Analysis 


Program 


%Ti-6Al-4V Yielding Analysis 
%Input the material properties 
EB 115.780: 
n= 9521: 

Dues] ie 

Sy=925e6; 
Sena 2 2 5 
Sres=-0.85* Sy 

saxıal EA 

YWhere E:Young’s Modulus, nu:Poisson’s Ratio, R:Radius of fretting pads. 
7 Sv: Yield Strength. Sres:Residual Stress, Saxial: Applied bulk stress 
Frage a 

mu 05: 

¡OA 

PSU 

Ry=0; 


% Increase the contact load, P. for each mu and Saxial until the effective 


f 
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% stress exceeds Sy. 

while Ry<1 
IR E R) (1/3): 

Ba SE, > pa 2) 
Eh: 

BZ mu Bie 

"Input Hamilton's Equations 
y=0; 

z(1)=0.0000001*a; 


Sl)= 1.092 


tor j=l] 
forse 
ات‎ 


A-(r^2)4(2()^2)-(a72) 
S—((A^2)-(1*(a^2)*(2()) ^2))) ^ (1/2) 

MENSTA OM 

N\=((S-A)/2) (1/2); 

phi=atan(a/M): 

G-(M^2)-(N^2) - (z(]) * M)-(a* N): 

HEN AN， 

%Stresses due to normal loads 

Sxn=(Po/a) (nuz) PRIT 2D) Oa EEN 
((1-2*mu)/3)* (N*S-E2* A*N--a^3)-nu* M*z(J)*a)-N lx(ı) لل‎ O AON 
Syn=(Po/a)*((1+nu)*z(j)*pli+(1/r"2)*((x(i) *2/r72)*((1-nu) * N*z(j) *2-((1-2* nu) /3) 
*(N*S4-2* A* N-Fa^3)-nu* M*z(j)*a)-N* (2*nu*x(1) ^2)-(0/8))); 
Szn=B1*(-N+(a*z(3)*M)/S); 

TAzxn-BI* -z(j)*(x()*N/S-(x()*z(j)* H/(G^2--H^2)))); 

%Stresses due to tangential loads 

Sxt=B2*(-x(i)* (nu/4+1)*phi+(a*x(i)*M /r74)* ((3/2-2*x(i) *2/r72)* (S*nu-2* A *nu+ 
Z(j)*2)+(x(G)*2*z(j)72/S)+(7* nu *r* 2 /4)-(2* nu*x (i) *2)491*2)4+(x(i)*2(j)*¥N/r74)* 


R] 





2 2/72)” E 2*nu)-A/3*(1-2*nu)-.5* (z()) 22-3*a^2))- 
(a^2*x(1) ^2/S)-(nu*r^2/4)-(r*r^2/4)) -(4*a^3*x(1)*z()/ (3*r^4))* 

me O 

Sytz B2*((-3*nu*x(1)*phi/4)-- (a*x(1)*M/r^4)*((1/2-2*0/r^2)* ($*nu-2* A*nuaz(j) ^2 
r^2*nu)2-(0*2()) ^2/S)-(3*nu*r^2/3)) -(x(1)*zG)*N/r^4)*((1/2-2*0/r^2)* 
(-8/6* (1-2*nu)-A/3* (1-2*nu)-.5* (z(]) ^2--3*a^2))2-(a^2*0/S)-(3*nu*r^2/4)-(1*r^2/4))4- 
pocos» ISSN Er OO 

1L DO ZU x0) iN /(2*672)* UTE (6^ 2:420) 22:527 2)7 5) 
TAzxt-B2*((3*z(j)*phi/2)--(a*z(j))* M/r^2)*(1--x(1)^2/r^2-x(1)^2/S) -N/r^2* 
(-.75*(S--2* A) -z(j) ^2-.15*a^2-.25*r^2-4-z(]) ^2/2* (.5-2*x(1)^2/r^2)); 

%Now combine the stresses 

Sxtot=Sxn+Sxt+5res+Saxial; 

Sytot=Svn+Svt+Sres; 

Sztot=Szn+Szt; 

TAzxtot=TAzxn+TAzxt; 

% Determine J2 and von Mises effective stress 
J2=1/6*((Sxtot-Sytot)”2+(Sytot-Sztot)”2+(Sztot-Sxtot)”2)+TAzxtot” 2; 

EB a 

Sen, 53 SEIS: 

x(i+1)=x(i)+.02*a; 

end 

z(j4-1)—z(]) 4-.025*a; 

end 

[Sem.i]=max(Se); 

GER 

xx]; 

2-2 

asa: 

Se=Se; 


%Now give x.z location of max stress 


RD 





[Smm.zm]=max(Sem); 

zuım=z(zm); 

zma=zınn/a; 

xm=i(zm); 

SONUS‏ اه 

xma=xmm/a; 

Now find the stresses at x=a, and x=-a at the surface 
E ERI 

Eo TPC pra 2) 

oy a. 

E Smu BI; 

po eu. 

p 

Sxtna-B2*(-pi*x2/2* ((nu/4)4-1)); 
syma=B2 e nu iS): 
Bass SD cm ao 

Sima Bl (124 ca 
Sxtona=Sxtna+Sxnna+Sres+Saxlal:; 
Sytona=Sytna+Synna+Sres; 
J2naz1/6* ((Sxtona-Svtona)^2--(Sytona-0) ^ 24- (0-Sxtona)^2) 4-0^2; 
H?nasgj2na .5: 

Sena=37.5*H2na; 


“%Find stresses at x=a 


Sxtpa=B2*(-pi*x2/2*((nu/4)+1)); 
Spa bod eU لل ا‎ S 
Spa 0015251) a) on. 
Synpa=B1*((1-2*nu)*(-a/3)); 


SxtoPa=Sxtpa+Sxnpa+Sres+Saxiıal; 


RA 





SytoPa=Sytpa+Synpa+Sres; 
J2pa=1 /6* ((SxtoPa-SytoPa) *2+(SytoPa-0)*2+(0-SxtoPa) *2)+072: 
H2pazJ2pa ^ .5; 

Depa.) Epa. 

%Now determine the maximum eflective stress 
Smmm=[Smm,Sena,Sepa] 
[Smmax.xmmax|=max(Smmm); 
Semax=5mmax; 


Ry—Semax/ Sy; 


=P +2: 
end 
ll 


P2{kJj)=P1; 

mu=mu-+.1: 

end 

Saxial=Saxial+(.1*Sy) 

end 

Emus 

pospo 

%Now put P terms over yield strength. 
[arm 13 

longer 

au Diva ka): 

Pori (3 E22 pita 2 
Bauen) Lollın) sy. 

end 


end 


Ryy=Ryy 


RA 










H 


الل 


f d 
NU 


seais O u fim D 
514 

Hu se | 

THETTA 


Appendix B 


Depth of Yielding Analysis 


%Ti-6Al-4V Depth of Yielding Analysis 
"Input the material properties 
5-1159 
Dur 
Kelle: 

Sv=925e6; 

EE 
Sres=-.89*Sy 
mu=0.3: 


sux Sy 


lore 
RETO: 
Ork? 


Se i le SE 

Eois pine, 
Io 

B2>= mu Bl. 

Z Input Hamilton's Equations 
y=0; 

z(1)=0.0000001*a; 





$ B 
INE. i n 
rw 






e EL 


TU 


Upon ves LL TN 


Se UO a, 

for JR 2 1 

ore حا‎ 1 !! 

م 

A-(r^2)4 (2()^2)-(a^2); 

S—((A^2)-(4* (a^2)*(26)^2))) (1/2); 

M=((S+A)/2)” (1/2); 

RAT a2: 

phi-atan(a/M); 

GES8N102)-(N^2)--(2()  M9-(a* N); 

H=(2*M*N)+(a*M)+(z(j)*N); 

STRESSES DUE TO NORMAL LOADS 
Sxn-z(Po/a)*((1-2-nu)*z(j)*phi--(1/r^2)*((0-x(i) ^2/r^2)* ((1-nu)* N*z(]) ^2- 
((1-2*nu)/3)* (N*S-2* A*N -a^3)-nu* M*z(])*a)-N* (x(1)^2)-(M*x(1) ^2*z(j)*a/8)): 
Sm oa) (lean) ai pice ra AN 2s 
((1-2*nu) /3)* (N*S+2* A*N+a73)-nu* M*z(j)*a)-N* (2* nu*x(i) ~2)-(0/5S))); 
SB Na ZO) MS 
MB ANS A AA 

“Stresses due to tangential loads 

SB Saa depues EE 
E A A ار‎ no 
SMN LL NL 

5*z(j) 2+3*a°2))+la°2*x(i)"2/S)-(nu*r"2/4)-(7’r2/4))+ 

网 2er 

EC EOC TU XU PHVA SUI MT ye po PT 

(S*nu-2* A*nu+z(j) "2+r"2*nu)+(0*z(j) "2/S)+(3*nu*r"2/4))+(xG)*z()N/r4)* 
(( 172-207 2)°(-5/6° (1-2mu)-A/ 3 12m. 32) ars: 
Stat LT AEE SS A ADO ele en) 
St =B2* ((a(j)*x(i)*N/(2*17-2)*(1-(r724-2(j)-2+a-2)/S))); 

TAzxt=B2* ((3*z(j)*phi/2)+(a*z(j)*M/r72)*(1+x(i) *2/r°*2-x(i) ©2/S)+N /r72* 


RA 





eget) 7) 2-.75*a 2.25% 57 2+2(j)°2/2*(.5-2* xl) "2/1 2))); 
%Now combine the stresses 

Sxtot=Sxn+Sxt+Sres+Saxial; 

Sytot=Syn+Syt+Sres: 

SA O= SZN OZ: 

TAzxtot=TAzxn+TAzxt: 

% Determine J2 and von Mises effective stress 

J2—1/6* ((Sxtot-Sytot) ^2-4- (Svtot-Sztot) ^ 2-- (Sztot-Sxtot) ^2)-- TAzxtot ^2; 
[7-12 «5: 

Set )=3 oH>: 

x(i+1)=x(1)+.02*a:; 

end 

z(j+1)=z(j)+.02* a; 

end 


[Sem.i=max(Se); 


a 
xlexl: 
21-77 


I la 

Se=Se; 

%Now give x.z location of max stress 
[Smm.zm]=max(Sem); 

AZ (ZE 

zma=zmmY/a; 

Ae 

AE 

xma=xmmY/a; 

“Now determine when von Mises effective stress equals yield stress 
Ra 


211192-42 


7 






"i: 
it. Va anilioob 9 p Ag 
AMO FSAI AA 


D 
i MË aa 


TIC 


١ 0 


> 


4 we 


alma a: 

x1(1)—-1.05*a; 

while Ry>1 

loro q^ pro 

a-(K*P*R)^(1/3); 

Po—(3*P)/ (2*pi*a^2); 

Bi 10,2 

Bon BE 

EXE 

Dal) 2), 

S=((A72)4+(4*(a72)*(2172)))" (1/2): 

M=((S+A)/2)" (1/2); 

N=((S-)/2)" (1/2); 

phi=atan(a/M): 

03125727131265: 

H=(2*M*N)+{a*M)+(z1*N); 

%Stresses due to normal loads 

Sn (Boa) ı (na) zZ pr 2er en Nl = 
((1-2*nu)/3)* (N*S--2* A*N--a^3)-nu* M*z1*a)-N*x1(1) ^2-(M*x1(i)^2*z1*a/8))); 
Sn IOI Inu ZI pn A Sa eN OS 
((1-2*nu)/3)* (N*S--2* A*N--a^3)-nu* M*z1*a)-N* (2*nu*x1(1) ^2)-(0/5))): 
sao ENA a ZMS): 

Ii ZI DENS ZH eS 

%Stresses due to tangential loads 

oxt- B2 cx Sina 4-1) *phissra ssp SNNT oS UL ل‎ 
(S*nu-2* A*nu4-z1^2)-(x1(1)^2*z1^2/S) 4-(7*nu*r^2/4)-(2*nu*x1(1)^2)-r^2) - 
Er EE Eege 
KE rte SIU) SOUT TT DISA a o ee 
Sy ein 


Syt=B? EE TEE 


RR 





r^2*mi)4-(0*z1^2/5)-(3*nu*r^2/4)) 4-(x1(1)*z1* N/r^ 4)*((1/2-2*0/r^2)*(-8/6* 
(1-2*g9-A/3* (1-2*nu)-.5* (z1L^24-3*a^2)) - (a^ 2*0/S)-(3*nu*r^2/4)-(1*r^2/4))4- 
eel) 71) (8514) (1 /2-2*0/r 2)*(1-2* nu); 

Sacer ZI EN (2*r^2)* CIe(r* 2-Ez1^2-Fa^2)/ 9) 

za = p250053*21 pn 2)--(a*zi*ND/r^2)* (I -exH(0)^2/r^2-x10)^2/8)-EN7r^2* 
i05 152 )5EZ132-50^5202-.280"p^2-E21^2/2*08:2*x 0) ^ 2/0 ^2) 34 

Now combine the stresses 

Sxtot=Sxn+Sxt+Sres+Saxial; 

Svtot—5vn4-Svt--Sres; 

Sztot=Szn+Szt; 

TAzxtot=TAzxn+TAzxt; 

^ Determine J2 and von Mises equivalent stress 

J221/6* ((Sxtot-Sytot)^2--(Sytot-Sztot) ^24- (Sztot-Sxtot) ^2)-- TAzxtot ^2; 
Bos 

Sola —395"H7 

x1(i+1)=x1(1)+.02%*a; 

end 

Seml=max(Sel); 

Ry=Seml/Sy 

a‏ ا 

end 

z11(k,3))=(z1/a)-.01; 

P=P+-0; 

end 

Saxial=Saxial+.1*Sy 

end 

Zll=zll 


Tu ra 
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Appendix C 


Maximum Principal Stress Analysis 


Ve Ti-6Al-4 V Maximum Principal Stress Analysis 
Input the material properties 
E=115.7e9; 
Du 00], 
] 273 
Sy=925e6: 
am oU 
N 3 e al: 
Sres=- 853 Sy; 


Baal | ov. 


(ee 
rn 10: 
iore — 


Exe 

a=(K PRI (I/33): 

Dor (pina: 

BIE Po a 

P= BI 

"Input Hamilton's Equations 


v=0; 


an 







ite P 
E 
o A 
p^ mm 
i48 


af 
ib A. i 











001 An rdi و‎ 


IN 





z=0.0000001*a; 

x(1)2-1.05*a; 

while Slmax>0 

%Slmax is the maximum principal stress 

[ord M 

rs 

ar 2122,02); 

S-((A^2)4-4*(a^2)*(2^2)))* (1/2); 

M=((S+A)/2)* (1/2); 

er ly: 

phi=atan(a/M); 

G-(M^2)-(N^2)-(z* M)-(a*N); 

H2(2*M*N)--(a* M) 4-(z* N); 

STRESSES DUE TO NORMAL LOADS 

an- roa) m stole ms» qc miae 
Eeer E a So usn 
n= EE E zie 

WM SAO ISLAN a So NE a ال ل‎ 
Szn=B1*(-N+(a*z*M)/S); 
EEN 

%Stresses due to tangential loads 

SLB CSU) ODI I DEA SI MT E O e a 
(S*nu-2* A*nu4z ESO) 22 SET RUT ETO SI) OTT 
ا‎ OT E Ta TaD 
OFT a la AS mar ee 2a 3 
Zio me (8/22 027722 m, 
Svt—B2*((-3*nu*x(1)*phi/4)-(a*x(0)*M/r^3)*(072-2*0/r^2)* 

(S*nu-2* A*nu4-z^2-Er^2*nu)J-(0*z^2/S)--(3*nu*r^2/4))--(x(3)*z*N/r^4)* 
((1/2-2*0/r^2)*(-S/6* (1-2*nu)-A /3*( I-2*m)-.5*(z^2-E3*a^2)) 4E 
(37220/5), mer A es 


Q1 





MN 2) (1-2*nu)); 

SZ; INE (Cz x(i)* \ /(2*r72)*(1-(r724272+a~2)/S))): 
gBuzscposa5z*pn/2)5(a*2* NI/r^2)* (1-Ex(1) ^2/r^2-x(1)^ 2/S]-F 
N/r^2*(-.15* (S--2* A)--z^2-.15*a^2-.25*r^2-pz^2/2* (.5-2*x(1)^2/r^2))); 
%Now combine tlie stresses 

Sutotsoxi-oxteb»les-oaxial: 

Sytot=Syn+Syt+Sres; 

Sztot=Szn+Szt: 

TAzxtot=TAzxn+TAzxt; 

ARO => NO 

SO V (1) =5Zt0L: 

TAzxtotv(ij=TAzxtot; 

% Determine the max principal stress for each point 
cxz=(Sxtot+Sztot)/2; 

rxzz((Sxtot-cxz) ^2-- TÀzxtot ^2) ^.5; 

SZ) cxz Lıxz: 

x(i+1)=x(i)+.02*a; 

end 

"Find the maximum principal stress of the vector of principal stresses 
SINE mas IZ je 

Exmaxesniaxi»xtoby 

Szmax=max(Sztotv); 

TAzxmax=max(TAzxtotv): 

Ma => 52m 

z=2+.01*a; 

zaa=z/a 

end 

22 cuz 0122: 

a=a; 


I= TOC 


P=P+40 


end 


Baile Saxıala.175Y 


end 
7-2 
Zu 


Inu-nm 


02 





Appendix D 
Hamilton’s Expressions 


D.1 Normal Load Stresses 





3P = A Se a ical 
Doc UD Ad ae 2 
Ita 7 _N (2? D 2u y?) 5 Maza 
2 2 $ 
3P IMs EE 
Qe TS / [K ) B ) | 
” Asa r2 


—N (y? + 2vr?) — nz E 








"m 3P |-N n | 


Ira? S 


At y = 0, the onlv non-zero shear stress Is 7;,. 


"XE zN xzH 
Ie | “| S Gg 


where 

€ € 一 a € / € € 
SE 
A SANZ oS en 
M= (>=); N= (>52); ó—tan (5) and 


OA 


i i then = (Ph ac Abg ف‎ ey.) AA Pv te 1 
1 glk  )اهدب اء‎ 


| i ۱ 
e i Wi / A 
f e 9 "te € i | 1 
| 0 A fi u 





pa 


G= Ai? — N24 2zAJ/ —aN: H=2NIN+aAl+2zN 
On the axis where r = 0. the formulae above are indeterminate and must be re-worked 


Hamilton performed tlis and the expressions are now given by: 


3P = E 
Oy = Oy = amas (1 +r) (z tan (a/ z) 一 a) 十 一 一 一 





and Try = Ty; 


Toa 


D.2 Tangential Load Stresses 




















—r (E - 1) Ó 十 EA HE 一 27) (Sv — 2Àv + 25) 一 — + au em r2} 
3Q D 942 LES E 272 
| FF) Emule 
DE 2 e 
He) 
2,2 
o + 一 一 iG = 2 | [(Sv — 2Àv + Do) kl “= svr” \ 
30 ; 


(1 —2v) - $1 = 2») - 1 (2? + 302)] + SE Sc 


4 
Aadzz [1 20€ e 
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Again the only non-zero shear stress at y — 0 is 7,4. 


* Note the equation for c; in the original paper bv Hamilton contained an error. The paper 


contained the number 3 in the expression marked above where it should have contained an $. 


The units do not work if this change 1s. not made. 
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